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Abstract 


The  design  and  generation  of  modified  involute  helical  gears  that  have  a  localized  and 
stable  bearing  contact,  and  reduced  noise  and  vibration  characteristics  are  described.  The 
localization  of  the  bearing  contact  is  achieved  by  the  mismatch  of  the  two  generating  surfaces 
that  are  used  for  generation  of  the  pinion  and  the  gear.  The  reduction  of  noise  and  vibration 
will  be  achieved  by  application  of  a  parabolic  function  of  transmission  errors  that  is  able  to 
absorb  the  almost  linear  function  of  transmission  errors  caused  by  gear  misalignment.  The 
meshing  and  contact  of  misaligned  gear  drives  can  be  analyzed  by  application  of  computer 
programs  that  have  been  developed.  The  computations  confirmed  the  effectiveness  of  the 
proposed  modification  of  the  gear  geometry.  A  numerical  example  that  illustrates  the  devel¬ 
oped  theory  is  provided. 


Nomenclature 


a  Parabola  parameter  (fig.  2(a)) 

Cc  Modification  coefficient  of  pinion  rack-cutter  surface  (fig.  5) 
b  Slope  of  linear  function  (fig.  1(b)) 

Epc  The  shortest  distance  between  the  pinion  and  rack-cutter  Sc  (fig.  6(b)) 
Epg  The  shortest  distance  between  the  pinion-gear  axes  (fig.  7) 

Egt  The  shortest  distance  between  the  gear  and  rack-cutter  St  (fig.  6(a)) 

77121  Gear  ratio 

Mtj  Coordinate  transformation  matrix  (from  Sj  to  5,) 
n.  Unit  normal  vector  to  rack-cutter  surface  Sr  (r  =  c,t) 

ny  ^  Unit  normal  vector  to  surface  S,-  represented  in  coordinate 

system  Sf  (i  =  p,sr) 

Ni  Number  of  teeth  of  the  pinion  (i  =  l,p)  or  the  gear  (i  =  2,g) 

Nr  Normal  vector  to  rack-cutter  surface  Sr  {r  =  c,t) 

Pn  Circular  pitch  in  normal  section  (fig.  3) 

P„  Diametral  pitch  in  normal  section 

r, -  Radius  of  the  pitch  circle  of  the  pinion  (or  gear)  (i  =  p,p) 

rt,r*  Position  vector  of  surface  S,- 

Ty  ^  Position  vector  of  surface  S,-  represented  in  coordinate  system  Sj 

s, Sr  Displacement  of  rack-cutter  Sr  (r  =  c,i)  (fig.  6) 

Si  Coordinate  system  i 

Ui,0i  Surface  parameters  of  S,-  {i  =  p,g) 

Ur, It  Surface  parameters  of  Sr  (r  =  c,  i) 

v(’’)  Velocity  of  rack-cutter  surface  point  (r  =  c,t) 

Relative  velocity  of  surface  S,-  point  with  respect  to  surface  Sj  point 
ao  Normal  pressure  angle  (fig.  3) 

/3(,  Helix  angle  on  the  pinion  (gear)  pitch  cylinder  (figs.  3  and  4) 

6  Elastic  approach  of  pinion  and  gear  tooth  surfaces 
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A.E  Change  of  center  distance 

AAo  Change  of  pinion  lead  angle  on  the  pitch  cylinder 

Aq  Displacement  of  contact  point  caused  by  misalignment 

A'yx  Misalignment  angle  formed  by  crossed  gear  axes  (fig.  8(a)) 

A7j,  Misalignment  angle  formed  by  intersected  gear  axes  (fig.  8(b)) 

A<^2AV’2  Transmission  error  (fig.  2) 

Vector  of  the  angle  of  compensating  turn  of  gear  2 
Lead  angle  on  pinion  pitch  cylinder 
Pinion  (i  =  p)  and  gear  (i  =  g)  tooth  surfaces 
Raxic-cutter  surfaces  (r  =  c,f) 

Rotation  angle  of  gear  i  (i  =  1,2,  p,g)  (figs.  2  and  7) 

Rotation  angle  of  gear  being  in  mesh  with  the  rack-cutter  St  (fig.  6(a)) 
Rotation  angle  of  pinion  being  in  mesh  with  the  rack-cutter  S^  (fig.  6(b)) 
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1.  Introduction 


Conventional  helical  involute  gears  are  designed  for  transformation  of  rotation  between 
parallel  axes.  Theoretically,  the  gear  tooth  surfaces  are  in  line  tangency  at  every  instant, 
along  a  straight  line  that  is  a  tangent  to  the  helix  on  the  gear  base  cylinder.  However,  the 
line  contact  of  gear  tooth  surfaces  can  be  realized  only  for  an  ideal  gear  drive.  In  reality,  the 
crossing  of  axes  of  rotation  (instead  of  being  parallel)  and  errors  of  lead  angle  result  in  the 
so-called  edge  contact,  as  a  specific  instantaneous  point  contact  caused  by  curve-to-surface 
tangency.  Here,  the  contacting  curve  is  the  edge  of  the  tooth  surface  of  one  of  the  mating 
gears  and  the  contacting  surface  is  the  tooth  surface  of  the  other  one. 

Trying  to  avoid  the  edge  contact,  the  manufacturers  of  helical  gears  use  various  methods 
of  crowning  (deviation)  of  the  theoretical  gear  tooth  surfaces.  However,  the  applied  methods 
of  crowning  have  not  been  complemented  with  the  analysis  of  transmission  errors  caused 
by  misalignment.  Our  investigation  shows  that  improper  crowning  may  avoid  edge  contact 
but  cannot  avoid  the  appearance  of  transmission  errors  of  the  shape  shown  in  fig.  1.  The 
function  of  sudi  transmission  errors  is  piecewise,  almost  linear,  and  has  the  frequency  equal 
to  the  cycle  of  meshing  of  one  pair  of  teeth.  The  above  mentioned  transmission  errors  cause 
high  vibration  and  noise  and  therefore  such  transmission  errors  must  be  avoided.  This  can 
be  achieved  by  application  of  computer  numerically  controlled  (CNC)  machines  that  have 
opened  new  perspectives  for  generation  of  gear  tooth  surfaces  with  improved  topology. 

The  intent  of  this  paper  is  to  describe  a  modified  topology  of  low-noise  involute  helical 
gears  that  satisfies  the  following  requirements: 

(1)  The  noise  and  vibration  of  helical  gears  are  reduced  substantially  by  application  of  a 
predesigned  function  of  transmission  errors  of  a  parabolic  type  (fig.  2).  Such  a  function  can 
absorb  (see  below)  an  almost  linear  function  of  transmission  errors  shown  in  fig.  1. 

(2)  The  bearing  contact  is  localized.  Theoretically,  the  tooth  surfaces  are  in  tangency  at 
every  instant  at  a  point  instead  of  a  fine.  The  contact  of  gear  tooth  surfaces  at  every  instant 
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is  spread  over  an  elliptical  area  due  to  elcistic  deformation  of  gear  teeth.  The  dimension  of 
the  instantaneous  contact  ellipse  can  be  controlled  by  choosing  proper  design  parameters. 

(3)  The  proposed  gear  tooth  surfaces  can  be  generated  by  two  rack-cutters  designed  for 
generation  of  the  pinion  and  gear,  respectively.  A  nonlinear  transmission  function  in  the 
process  for  gear  generation  must  be  provided  and  this  can  be  accomplished  by  application  of 
the  CNC  machine.  A  linear  transmission  function  is  provided  in  the  process  for  the  pinion 
generation. 


2.  Interaction  of  Parabolic  and  Linear  Function  of  Transmission  Errors 


Ni 


The  ideal  gears  transform  rotation  with  constant  gear  ratio  m2i  =  and  the  ideal 

JV2 

transmission  function  is 


where  Ni  and  are  the  tooth  numbers  of  the  pinion  and  gear,  respectively. 

However,  the  crossing  of  gear  axes  (instead  of  being  parallel),  intersection  of  these  axes 
and  errors  of  lead  angle  cause  a  transmission  function  ^2(^1)  that  is  shown  in  fig.  1(a). 
Our  investigation  (see  sections  4-6)  shows  that  the  function  of  transmission  errors  caused  by 
above  mentioned  errors  of  misalignment  is  a  piecewise  almost  linear  function  of  transmission 
errors  A<f>2{<l>i)  with  the  frequency  of  a  cycle  of  meshing  for  one  pair  of  teeth  (fig.  1(b)). 

Here: 


(2) 


Transmission  errors  of  this  type  cause  a  discontinuity  of  the  transmission  function  and  a 
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big  jump  of  the  angular  velocity  of  the  driven  gear  at  transfer  points  (when  one  pair  of  teeth 
is  changed  to  another  one).  Therefore,  vibration  and  noise  become  inevitable. 

It  was  proven  [1,2,4]  that  a  predesigned  parabolic  function  of  transmission  errors  in¬ 
teracting  with  a  linear  function  will  become  a  parabolic  function  with  the  same  parabola 
parameter.  A  parabolic  function  of  transmission  errors  is  much  more  preferable  than  a  lin¬ 
ear  function  since  the  transmission  function  will  be  a  continuous  one,  the  jump  of  angular 
velocity  of  the  driven  gear  and  the  stroke  at  the  transfer  point  will  be  substantially  reduced. 

Fig.  2(a)  shows  the  sum  of  two  functions  of  transmission  errors 

A^2(^i)  —  ^^2  "b  ^^2  ^(^i)  ~  ~ 


The  first  one,  A<^2^^(^i),  is  caused  by  misalignment.  The  second  one,  A<f)^\<j>i),  is  a 
predesigned  parabolic  function  which  exists  even  if  misalignment  does  not  appear.  It  is  easy 
to  verify  that  equation  (3)  represents  in  the  new  coordinate  system  (At/’2,  ^i)  the  parabolic 
function  (fig.  2(b))  that  is  designated  as 

A7p2  =  — aV’i 


The  parabola  parameter  a  in  equations  (3)  and  (4)  is  the  same.  Axes  of  coordinate 
system  (AV>25  V’l)  and  (A<^2,  <i>i)  are  parallel  but  the  origins  are  different.  The  coordinate 
transformation  from  (A<^2,  <^1)  to  (Ai/>2,  ^1)  is  represented  with  the  following  equations 


b 

AV’2  =  -  — 


(5) 
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The  difference  between  functions  A<f>2{<i>i)  and  AV’2(V»i)  is  the  location  of  the  couple  of 
points  {A,B)  and  the  respective  points  (A*,  B*)  (fig.  2(a)).  The  symmetrical  location  of 
(A,B)  is  turned  into  the  asymmetrical  location  of  (A*,B*).  However,  the  interaction  of 
several  functions  AV»2(V’i)  determined  for  several  tooth  surfaces  being  in  mesh  may  provide 
a  symmetrical  parabolic  function  of  transmission  errors  as  shown  in  fig.  2(b).  This  can  be 
achieved  if  the  parabolic  function  Ahi^i)  will  be  predesigned  in  the  area  (fig.  2(a)) 


MB)  -  MA)  >  ^  +  2c 


(6) 


where  c  =  — .  Requirement  (6),  if  observed,  enables  to  provide  a  continuous  function 

2a 

AV’2(t(’i)  for  the  range  of  ^  where  Ni  is  the  pinion  tooth  number.  It  will  be  shown  below 
(see  sections  5  and  6)  that  functions  of  transmission  errors  caused  by  angular  errors  (such 
as  the  crossing  and  intersection  of  the  axes  of  rotation,  error  of  the  lead  angles)  are  indeed 
piecewise  linear  functions,  and  the  coefficient  b  can  be  determined  knowing  the  angular  error 
caused  by  misalignment  and  the  design  parameters  of  the  gear  drive. 


3.  Surfaces  of  Rack-Cutters 

The  imaginary  process  of  generation  of  conjugate  tooth  surfaces  is  based  on  application 
of  two  rack-cutters  that  are  provided  respectively  by  a  plane  S*  and  a  cylindrical  surface  Sc 
that  differs  slightly  from  plane  S(  (see  fig.  3).  The  rack-cutter  surfaces  St  and  Sc  are  rigidly 
connected  each  to  other  in  the  process  of  the  imaginary  generation,  and  they  are  in  tangency 
along  a  straight  line,  ObZh  (fig.  5).  This  line  and  the  parallel  axes  of  the  gears  form  angle 
M  that  is  equal  to  the  helix  angle  on  the  pinion  (gear)  pitch  cylinder.  The  normal  sections 
of  the  rack-cutters  are  shown  in  figs.  3  and  5.  Rack-cutter  surface  Sc  generates  the  pinion 
tooth  surface  Sp,  and  the  rack-cutter  surface  St  generates  the  gear  tooth  surface  Sp. 
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Gear  Rack-Cutters  Ej 


Using  figs.  3,  4  and  5,  we  represent  the  transformation  matrix  from  system  Sa  to  Sr 
(r  =  c,  t)  and  Sb  to  5a  as  follows 


Mr.= 


10  0  0 

0  cos^o  sin  00  0 

0  —  sin^o  cos  00  0 

0  0  0  1 


(7) 


ab  — 


cos  Go 

no 
0 
0 


sin  Oo  0  —dp  cos  Oo 
0  am— dcsinao 
0  1  0 

0  0  1 


(8) 


Mrb  = 


COS  Oo  —  sin  Uo  0  —dp  cos  oco 

sin  ao  cos  0o  cos  Oo  cos  0o  sin  0o  (om  —  dp  sin  a<,)  cos  0o 
—  sin  OLo  sin  0o  —  cos  a,  sin  0o  cos  0o  —{dm  ~  dp  sin  a^)  sin  0o 
0  0  0  1 


(9) 


Here  Oo  and  0o  are  the  normal  pressure  angle  and  the  helix  angle  of  the  rack-cutter; 
is  the  half  of  the  tooth  width  of  the  rack-cutter  on  middle  line  m-m  (fig.  3),  where 


TT 

TK 


(10) 
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and  P„  is  the  normal  diametral  pitch  of  the  rack-cutter,  is  the  distance  between  middle 
line  m  -  m  of  the  rack-cutter  and  the  origin  Oh  along  axis  xj,  fig.  5.  Parameter  dp  can  be 
controlled  to  adjust  the  location  of  the  contact  path  on  the  gear  tooth  surface. 

Surface  Et  of  the  gear  rack-cutter  is  a  plane  that  is  represented  in  Sb  as 

J  ut  0  It  I 


where  (ut,  h)  are  the  surface  parameters. 

Rack-cutter  surface  Et  is  represented  in  coordinate  system  St  by  the  matrix  equation 


rt(uf,  h)  =  Mt6rf  ^ 


(12) 


Equations  (9),  (11)  and  (12)  yield 


rt 


(ttt  —  dp)  cos  ag 

[(ut  -  dp)  sina^  +  a„,]  cos  +  h  sin  0g 
—  [(ut  —  dp)  sin  ag  -b  a,„]  sin  /?<>  -t-  h  cos  ^g 


(13) 


The  unit  normal  to  Et  is  represented  in  5t  by  equations 


Nt  M  ^  V 
”  INtl  ’  ‘  dlt  dut 


(14) 


that  yield 
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(15) 


nt  =  —  sin  Oo  cos  cos  —  cos  sin  0o 

Pinion  Rack-Cutter  Surface  Sc 

Rack-cutter  Ec  generates  the  pinion.  The  normal  section  of  rack-cutter  surface  Sc  (fig.  5) 
is  a  parabolic  curve.  We  remind  that  the  normal  section  of  rack-cutter  surface  S(  is  a  straight 
line  directed  along  axis  Xb  in  fig.  5.  The  pajabolic  curve  is  in  tangency  with  the  zj-axis  at 
point  N{Ob)‘  Rack-cutter  surfaces  Sc  and  S*  are  in  tangency  along  a  straight  line  that  is 
parallel  to  axes  and  Zb  and  passes  through  point  Ob  that  coincides  with  point  N.  The 
deviation  of  the  parabolic  curve  from  the  x^-axis  affects  the  dimensions  of  the  instantaneous 
contact  ellipse. 

Rack-cutter  surface  Sc  is  represented  in  Sb  as  follows 

rl‘'=[u,  -a,u?  kf  (16) 

where  Oc  is  coefficient  of  the  parabolic  normal  section,  and  (uc,  h)  are  the  surface  parameters 
of  Sc- 

Rack-cutter  surface  Sc  is  represented  in  coordinate  system  Sc  by  the  matrix  equation 

r.(t<„y  =  Mdr<‘l  (17) 

Equations  (9),  (16)  and  (17)  yield 
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(uc  —  dp)  cos  Go  +  Ucul  sin  ao 


Tc  = 


[(uc  —  dp)  sin  a,,  +  am]  cos  0o  ~  o,c^\  cos  oco  cos  8o  +  h  sin  /3o 


-[(uc  —  dp)  sin  ao  +  a„]  sin  fSo  +  aeu\  cos  Oo  sin  /3<,  +  h  cos  jdo  J 


(18) 


The  unit  normal  of  Sc  is  represented  as 


„  duo  ^  dh 

'  l^X  ^ 

^duo  dk 


(19) 


Equations  (18)  and  (19)  yield 


1 

(l+4a2u2)0.5 


sin  Oo  —  SccUc  cos  ao 


—(cos  ao  +  2acUc  sin  o:,,)  cos 


(cos  ao  +  2acUc  sin  ao)  sin  /3o 


(20) 


Using  equations  (13),  (15),  (18)  and  (20),  it  is  easy  to  verify  that  surfaces  Sc  and  St  are 
in  tangency  along  the  Zh  axis  when  Ug  =  ut  =  0. 

4.  Pinion  and  Gear  Surfaces  Generated  by  Rack-Cutters 

In  the  process  for  generation  the  two  rigidly  connected  rack-cutters  perform  translational 
motion  while  the  pinion  and  the  gear  perform  rotational  motions  as  shown  in  fig.  6.  To 
provide  a  predesigned  parabolic  function  of  transmission  errors  for  each  cycle  of  meshing,  it 
is  necessary  to  observe  certain  relations  between  the  motions  of  the  rack-cutters  and  gears, 
respectively. 
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The  angle  ^pc  of  pinion  rotation  and  the  displacement  Sc  of  rack-cutter  Sc  are  related  by 
the  following  linear  function 


'fl’pc  = 


Sc 


(21) 


Here:  Vp  is  the  radius  of  the  pinion  pitch  cylinder. 

The  angle  ^gt  of  gear  rotation  and  the  displacement  $t  of  rack-cutter  St  are  related  as 
follows 


=  #{r)  - 

Here:  Np  and  Ng  are  the  tooth  numbers  of  the  pinion  and  gear,  respectively,  and  is 
the  initial  position  angle  of  the  gear  for  the  modification  gear  rotation. 

Equation  of  Meshing  between  Rack-Cutter  Sc  and  Pinion  Sp 

The  equation  of  meshing  between  rack-cutter  Sc  and  the  pinion  tooth  surface  Sp  is 
represented  as 


/K,  h,  v-,c)  =  Nfl  •  vW  =  0 


(23) 


where  V’pc  is  the  angle  of  rotation  of  the  pinion  in  the  process  for  generation.  The  normal 
to  Sc  in  Sc  can  be  obtained  by  equation  (20),  and  the  relative  velocity  of  the  pinion 
with  respect  to  Sc  may  be  represented  as 
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(24) 


X  (Rp  +  Fc)  -  (  0  Tp  0  )^ 


dt 


Here:  Rp  =  (OpOe)c  =  (^p  r^i^pc  0)^,  Tp  =  Ej,c  is  the  radius  of  the  pitch  cylinder  of 
the  pinion  (fig.  6),  0  0  1  )^. 

Substitution  of  equations  (19)  and  (24)  into  (23),  yields  the  following  equation  of  meshing 
between  Sp  and  Sc 


/(tic,  /c,  V’pc)  =  h  sin  00  +  Tplkpc  +  Om  COS  0o  + 


[(uc  -  dp)  +  2alul]  cos  0o 
sin  Oo  —  2acUc  cos  ao 


=  0 


(25) 


Surface  of  Pinion  Sp 

In  the  process  of  generation  of  pinion  surface,  rack-cutter  Sc  performs  uniform  translation 
and  the  workpiece  performs  uniform  rotation  (fig.  6(b)).  The  transformation  matrix  from 
system  Sc  to  Sp  can  be  represented  as 


Mpc  = 


cos^pc  sint/>pc  0  rpCosV>pc +  rpV’pcsinV’pc 
-  sin  V»pc  cos  xppc  0  -fp  sin  ^pc  4-  Tpij}pc  cos  ^pc 
0  0  1  0 

0  0  0  1 


(26) 


Pinion  surface  Sp  in  system  Sp  is  represented  as 


^p(nc5  '^pc}  hdpcFc 


Ic  =  -{-7^ — cos  0o  +  rp^pc  +  Cm  COS  0o} /  sin  00 
sm  Oo  —  ZOcUc  cos  Qo 


(27) 
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Substituting  equation  (26)  into  (27),  we  obtain  equation  of  Sp  as 

i-p(«c,V’pc)  =  i*p(«p,^p)  (28) 

Equation  of  Meshing  between  Rack-Cutter  St  and  Gear  S^ 

The  equation  of  meshing  between  rack-cutter  Sf  and  the  gear  tooth  surface  Sj  is  repre¬ 
sented  as 

/(«.,  I.,  V-,.)  =  n!‘>-v!’"  =  0  (29) 

where  xpgt  is  the  angle  of  rotation  of  the  gear  in  the  process  for  generation.  The  unit  normal 
n!‘^  to  St  in  St  is  represented  by  equation  (15),  and  Vt^*^  is  the  relative  velocity  of  the  gear 
with  respect  to  rack-cutter  St. 

We  recall  that  the  rack-cutter  St  performs  translation  with  constant  velocity,  but  the 
gear  performs  rotation  with  variable  angular  velocity  that  is  represented  as  (see  equation 
(22)) 

(30) 

The  relative  velocity  is  represented  as 

vW  =  X  (R,  +  r.)  -  (  0  0  (31) 
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where 


R3  =  (  -rg 


iV„ 


■r.rl^pc  0  f 


(32) 


and  rg  =  Egt  is  the  radius  of  pitch  cylinder  of  the  gear  (fig.  6). 

Substitution  of  equations  (15),  (30),  (31)  and  (32)  into  (29)  yields  the  following  equation 
of  meshing  between  Sg  and  S* 


Nr, 


f{uu  lu  ^gt)  =  (“t  -  dp)  COS  /9o  +  sin  ao{lt  sin  /?<>  +  j^rgippc  +  am  cos  /?„) 
2aNgrg{t{;pc  -  0^°^)  ^ 

— - ;  Y  :  : - -7^  cos  Go  cos 

iVp  -  2aNg{ippc  - 


(33) 


+ 


where 


>/>ji  =  -  “(V’s.  - 


(34) 


Surface  of  Gear  Eg 

It  must  be  remembered  that  the  gear  with  the  tooth  surface  Eg  performs  rotation  about 
its  axis  with  varied  angular  velocity  while  rack-cutter  Ej  performs  uniform  translation  (fig. 
6(a)).  The  transformation  matrix  from  system  St  to  Sg  can  be  represented  as 
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(35) 


N 

-  cos  xl^gt  sin  i>gt  0  rgCosxl}gt  +  rgj^tl)pcSinipgt 

—  sin  x/^gt  —  cos  xj^gt  0  Vg  sin  Xpgt  —  cos  xj^gt 


0 

0 


0 

0 


1 

0 


0 

1 


Gear  tooth  surface  Eg  in  system  Sg  can  be  represented  as 


It  -  -{(ut  -  > 


I  2aJV'grj(^6pe  cot^l 


(36) 


The  derivation  of  equation  (36)  is  based  on  transformation  of  equations  (33)  and  (34). 
Equations  (35)  and  (36)  enable  to  represent  the  gear  tooth  surface  in  two-parameter  form 
as  follows 


'Pgi'^t^'^gt)  —  “P  gi}^gi^g) 


(37) 


5.  Computerized  Simulation  of  Meshing  and  Contact  of  Pinion-Gear  Tooth 
Surfaces 

We  consider  that  the  surfaces  of  the  pinion  and  the  gear  generated  by  worms  Eu,  and  E;, 
are  represented  in  coordinate  systems  Sp  and  Sg,  respectively.  The  fixed  coordinate  system  5/ 
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is  rigidly  connected  to  the  housing  of  the  gear  drive  (figs.  7  and  8).  The  movable  coordinate 
systems  Sp  and  Sg  are  rigidly  connected  to  the  pinion  and  the  gear,  respectively.  An  auxihary 
coordinate  system  Sh  is  appUed  for  simulation  of  meshing  when  the  gear  axis  is  crossed  or 
intersected  with  the  pinion  axis  instead  of  being  parallel,  and  when  the  shortest  distance 
between  the  pinion  and  gear  axes  is  changed.  The  errors  of  misalignment  are  referred  to  the 
gear.  The  misalignment  angle  A-y  is  decomposed  into  two  components,  Afx  and  Afy  that 
represent  the  crossing  angle  and  the  intersection  angle,  respectively.  The  pinion  performs 
rotational  motion  about  the  2:/-axis.  The  axis  of  gear  rotation  is  Zh.  The  shortest  distance 
between  the  axes  of  rotation  is  designated  as  Epg. 

The  rotation  matrices  from  system  Sh  to  Sf  for  crossed  and  intersecting  angles  are  rep¬ 
resented  in  the  followings  (fig.  8) 


Lfh 


-10  0 

0  —  cos  A%  sin  A'jx 

0  sin  A7x  cos  A7x 


(38) 


Lfk 


—  cos  A7y  0  —  sin  A-ify 

0-10 

—  sin  A^y  0  cos  A-yy 


(39) 


We  represent  the  pinion  and  gear  tooth  surfaces,  Ep  and  Sj,  and  their  unit  normals  in 
coordinate  system  Sf.  The  conditions  of  continuous  tangency  of  surfaces  Sp  and  Sj  are 
represented  by  the  following  equations  [1,2]. 


rf{up,6p,<f>p) 


(40) 
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njrVp.  0p,  (l>p)  =  n^f\ug,  Og,  4>g) 


(41) 


Vector  equation  (41)  provides  only  two  independent  equations  since  =  |njf^|  =  1. 
The  total  number  of  independent  equations  provided  by  (40)  and  (41)  is  five  that  relate  six 
parameters 


/j(Up,  0p,  <j>p,  Ug,  9g,  (f>g)  -  0  (z  -  1 , 2,  .  .  .  5) 


(42) 


The  continuous  solution  of  the  system  of  nonlinear  equations  (42)  is  based  on  the  following 
procedure: 

(1)  Using  an  initial  guess,  we  determine  a  set  of  parameters  that  satisfy  equation  system 
(42).  Thus 


p(O)=(„(0)^(,(0)^„(0)_M0)_J0)) 


(43) 


(2)  One  of  the  variable  parameters,  say  <j>p,  is  chosen  as  the  input  one,  and  is  supposed 
that  the  Jacobian 


D(h,h,h,Si,h) 

D{Up^  Ug^9g^<l>g) 


(44) 


differs  from  zero.  The  derivatives  in  the  Jacobian  are  taken  at  point 

(3)  Then,  equation  system  (42)  can  be  solved  in  the  neighborhood  of  by  functions 
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“p(^p)i  ^p(^p)>  “s(<^p)>  Wp) 


(45) 


(4)  Vector  function  rp(up,0p)  that  determines  the  pinion  surface  Sp  and  functions  Up(<^p), 
dp{<f>p)  enable  to  determine  the  path  of  contact  on  Ep. 

(5)  Similarly,  we  can  obtain  the  path  of  contact  on  the  gear  surface  Sj  using  vector 

function  Tg{ug,0g)  and  functions  Ug{<f>p),  9g{4>g). 

(6)  The  paths  of  contact  on  pinion  and  gear  tooth  surfaces  slightly  deviate  from  helices 
in  the  case  of  an  aligned  gear  drive.  The  line  of  action  for  an  aligned  gear  drive  (the  set  of 
contact  points  in  Sj)  slightly  deviates  from  a  straight  line  that  is  parallel  to  the  gear  axes. 

(7)  The  transmission  function  <l>g{<f>p)  deviates  from  the  ideal  transnussion  function,  and 
the  function  of  transmission  errors  coincides  with  the  predesigned  parabolic  function. 

(8)  The  determination  of  dimensions  and  orientation  of  the  instantaneous  contact  ellipse 
needs  the  knowledge  of  the  principal  curvatures  and  directions  of  contacting  surfaces  and  the 
elastic  approach  of  surfaces.  This  problem  can  be  substantially  simplified  if  the  pinion-gear 
principal  curvatures  and  directions  are  expressed  in  terms  of  the  principal  curvatures  and 
directions  of  the  generating  surfaces  and  parameters  of  motion  [1,2]. 

6.  Numerical  Example 

The  method  developed  in  this  report  is  illustrated  with  the  example  discussed  below.  The 
design  parameters  of  the  pinion  and  gear  are  listed  in  Table  1.  The  numerical  simulation 
of  meshing  is  performed  for  an  aligned  and  misaligned  gear  drives  with  various  errors  of 
alignment  for  the  pinion  and  gear. 

Case  1.  Aligned  gear  drive 

Figure  9  shows  the  transmission  errors  for  the  aligned  gear  drive.  The  TCA  performed 
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Table  1:  Design  parameters  of  pinion  and  gear 


pinion 

gear 

tooth  number 

0 

CM 

II 

Ng  =  100 

1 

r,  r  1 

normal  diametral  pitch 

Pn=5  — 
in 

Fn  —  b  — 
in 

normal  pressure  angle 

ao  =  20" 

ao  =  20" 

helix  angle  on  pitch  cylinder 

^0  =  30" 

0 

0 

CO 

II 

tooth  length 

L=1.6  in. 

L=1.6  in. 

modification  coefficient 

Cc  =  0.0008 

a  =  0.0014 

elastic  approach 

S  =  0.007  mm. 

S  =  0.007  mm. 

confirms  that  the  predesigned  parabolic  function  of  transmission  errors  exists.  The  maxiraum 
transmission  error  is  approximately  8  axe  seconds.  Figure  10  shows  the  contact  pattern  and 
the  contact  path.  The  path  of  contact  on  the  tooth  surface  is  in  the  longitudinal  direction. 
The  major  axis  of  the  instantaneous  contact  ellipse  is  6  mm  for  the  assumed  elastic  approach 
of  the  surfaces  equal  to  0.007  mm. 

Case  2.  The  pinion-gear  rotation  axes  are  crossed 

Figures  11  and  12  show  the  transmission  errors  and  the  contact  pattern  for  the  case  when 
the  crossing  angle  A'Tctoss  is  4  arc  minutes.  The  maximum  transmission  error  is  8  arc  seconds 
and  contact  paths  are  shifted  up  and  down  on  the  gear  and  pinion  surfaces,  respectively. 
Figures  13  and  14  show  the  transmission  errors  and  contact  pattern  for  A-y cross  =  ~4  arc 

minutes. 

Case  3.  The  pinion-gear  rotation  axes  are  intersected 

Figures  15  and  16  show  the  transnaission  errors  for  the  misalignment  A7, nteraect  =  4 
arc  minutes.  Figures  17  and  18  show  the  transmission  errors  and  contact  pattern  for  the 
mentioned  above  error  of  alignment. 
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Case  4.  Influence  of  error  of  lead  angle 

Figures  19  and  20  show  the  transmission  errors  and  contact  pattern  when  the  error  of 
the  lead  angle  is  4  arc  minutes. 

Figures  21  and  22  show  the  transmission  errors  and  contact  pattern  when  the  error  of 
the  lead  angle  is  —4  arc  minutes. 

For  all  of  above  cases,  the  maximum  transmission  error  does  not  exceed  8  arc  second 
(with  very  small  deviations  of  this  value). 

7.  Conclusion 

From  the  analytical  study  presented  in  this  report  the  following  conclusions  can  be  drawn: 

(1)  The  interaction  of  a  parabolic  and  a  linear  functions  of  transmission  errors  has  been 
discussed  to  prove  the  possibility  to  absorb  almost,  linear  functions  of  transmission 
errors  caused  by  misalignment. 

(2)  Mismatched  surfaces  of  two  rack-cutters  for  generation  of  modified  involute  gears  have 
been  proposed. 

(3)  Generation  and  geometry  of  pinion-gear  modified  tooth  surfaces  have  been  determined. 

(4)  Computerized  simulation  of  meshing  and  contact  of  pinion-gear  tooth  surfaces  has  been 
developed. 

(5)  An  algorithm  for  determination  of  relations  between  the  curvatures  of  the  generating 
and  the  generated  surfaces  has  been  developed. 

(6)  An  algorithm  for  determination  of  the  contact  elhpse  has  been  developed. 

(7)  Directions  for  users  of  application  of  developed  computer  programs  for  the  design  of 
gears  with  the  modified  geometry  and  computerized  simulation  of  their  meshing  and 
contact  have  been  developed  (Appendix  C). 
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Parallel  to  gear  axes 


Fig.  4  Orientation  of  rack-cutters  with  respect  to  gear  axes 
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(a) 


Fig.  8  Coordinate  systems  and  Sf 
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Tpn^smssio^l  errors  J  second  l 


(NO  RSSEMBU  ERRORS} 


PINION  ROTRTION  ANGLE  (degree) 


Fig.  9  Transmission  errors  for  the  aligned  gear  drive 
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Fig.  10  Contact  paths  and  pattern  for  the  aligned  gear  drive 
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CROSSING  RNGIX  BETWEEN  AXES:  ^  MINUTES 


Fig.  11  Transmission  errors  for  the  alignment  error 
ATcross  =  4  arc  minutes 
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CROSSING  RNGLE  BETWEEN  AXES:  4  MINUTES 


Fig.  12  Contact  pattern  for  the  alignment  error 
ATcross  =  minutes 
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TRflNSniSSION  ERRORS  I  second) 


CROSSING  RNGLE  BETWEEN  AXES:  -i  MINUTES 


Fig.  13  Transmission  errors  for  the  alignment  error 
ATcross  =  -  4  arc  minutes 
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CROSSING  RNGLE  BETWEEN  RXES:  MINUTES 


RXIS  2  MM 


Fig.  14  Contact  pattern  for  the  alignment  error 
^ Across  =  “  4  arc  minutes 


TRANSMISSION  ERRORS  I  second) 


INTERSECTING  RNGLE  BETWEEN  fiXES;  4  MINUTES 


Fig.  15  Transmission  errors  for  misaligned  gear  drive  with 
ATi^tersecuing  =  ^  arc  minutes 
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INTERSECTING  ANGLE  BETWEEN  AXES:  MINUTES 


Fig.  17  Transmission  errors  for  the  misaligned 
gear  drive  with  ATi^tersectung  =  “  4  arc  minutes 
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TOOTH  HEIGHT  MM 


intersecting  RNGLE  between  AXES:  “4  MINUTES 


Fig.  18  Contact  pattern  for  the  misaligned 
gear  drive  with  A7ii,tersecttmg  =  "  4  arc  minutes 
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TRRNSniSSlON  ERRORS  I  second I 


ERROR  OF  LEAD  RNGLE: 


4  MINUTES 


Fig.  19  Transmission  errors  for  the  misaligned 
gear  drive  with  AXq  =  4  arc  minutes 
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ERROR  OF  LERD  RNGLE:  4  MINUTES 


AXIS  Z  MM 


Fig.  20  Contact  pattern  and  contact  path  for  misaligned 
gear  drive  with  AXq  =  "1  £^rc  minutes 
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TRANSMISSION  ERRORS  I  second I 


ERROR  OF  LERD  RNGLE;  “4  MINUTES 


Fig.  21  Transmission  errors  for  the  misaligned 
gear  drive  with  AXq  “  “  minutes 
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ERROR  OF  PRESSURE  RNGLE:  -I  MINUTES 


AXIS  Z  MM 


Contact  pattern  and  contact  paths  for  the  misaligned 
gear  drive  with  AXq  “  ”  minutes 


Appendix  A.  Relations  between  the  Curvatures  of  the  Generating  and  Generated 
Surfaces 

Direct  Relations  between  Principal  Curvatures  and  Directions  of  Mating  Sur¬ 
faces 

The  main  advantage  of  this  approach  (proposed  by  Litvin)  is  the  possibility  to  determine 
the  principal  curvatures  and  directions  of  the  generated  surface  in  terms  of  principal  curva¬ 
tures  and  directions  of  the  generating  tool  surface,  and  the  parameters  of  motion.  In  this 
case,  the  tool  surfaces  for  the  generation  of  the  gear  and  the  pinion  tooth  surfaces  are  rack 
cutters.  The  equations  developed  permit  a  simplified  computational  procedure. 

The  system  of  equations  that  relate  the  principal  curvatures  and  directions  of  the  gen¬ 
erating  and  generated  surfaces  is  as  follows.  Consider  that  unit  vectors  e/  and  represent 
the  principal  directions  on  the  tool  surface  2^  at  point  P  of  tangency  of  surfaces  Si  and  S2 
(fig.  Al).  The  principal  curvatures  on  the  mating  surfaces  /cy  and  kh  of  the  tool  are  given; 
the  parameters  of  motion  (see  below)  are  also  given. 

The  goal  is  to  determine  angle  cr  that  is  formed  by  unit  vectors  ey  and  e^,  and  principal 
curvatures  k,  and  k,.  (The  unit  vectors  e,  and  e,  represent  the  principal  directions  on 
surface  E2).  The  system  of  equations  for  determination  of  a,  Ks,  and  Kg  is  as  follows. 


tan  2(7  = 


_ 2615625 _ 

^25  ~  ^15  ~  (^/  “  '^/i)^33 


2615625  _  625  —  6j5  (/Cy  ^fe)t33 

t33  sin  2cr  fas  cos  2cr 


Here: 


=  «/  +  + 


615  +  625 

<33 


(Al) 

(A2) 

(A3) 
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bis  =  •  e/j)  - •  e/) 

625  =  (u;(^2)-e/)-Kfe(v(^^^-efc)  (^5) 

■  <33  =  -n  •  [(u;(^)  X  vj?)  -  x  vj^^)] 

(A6) 

4-(n  X  •  e/)^  - 

The  nomenclature  for  equations  (A4)  to  (A6)  is  described  as  follows: 

angular  velocity  of  the  generating  tool 
.  angular  velocity  of  the  generated  gear 
defined  as 

transfer  motion  velocity  of  the  generating  tool 
transfer  motion  velocity  of  the  generated  gear 
defined  as 

n  surface  unit  normal  vector 

The  equations  discussed  above  are  used  in  the  TCA  program  for  determination  of  the 
contact  ellipse  at  the  points  of  contact  path  of  the  modified  helical  gear  drive. 

Numerical  Example 

The  input  and  output  for  the  determination  of  the  principal  curvatures  of  the  pinion 
tooth  surface  are  shown  in  Tables  Al  and  A2.  The  input  and  output  for  the  determination 
of  the  principal  curvatures  of  the  gear  tooth  surface  are  shown  in  Tables  A3  and  A4. 
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Table  A1  Input  Data 


Description 

Symbol 

Values 

first  principal  curvature  of  tool 

Kf 

-0.0016  (1/mm) 

second  principal  curvature  of  tool 

first  principal  direction  of  tool 

e/ 

[0.0016  -  0.4999  0.8660]^ 

second  principal  direction  of  tool 

efc 

[0.9387  0.2992  0.1710]^ 

angular  velocity  of  tool  (1/sec) 

c,(i) 

[0.0  0.0  0.0]^ 

angular  velocity  of  pinion  (1/sec) 

u;(2) 

[0.0  0.0  1.0]^ 

transfer  velocity  of  tool 

(mm/sec) 

vg' 

[0.0  58.6588  0.0]^ 

transfer  velocity  of  pinion 

(mm/  sec) 

vg' 

[-0.1514  58.5951  O.Of 

surface  normal  of  tangent  point 

n 

[-0.3420  0.8138  0.4698]^ 

Table  A2  Output  Data 


Description  (for  pinion) 

Symbol 

Values 

first  principal  curvature 

-0.001543  (l/mm) 

second  principal  curvature 

K, 

0.03907  (l/mm) 

first  principal  direction 

e. 

[0.1751  -  0.4361  0.8827p 

second  principal  direction 

e. 

[0.9222  0.3865  0.0080]^ 

47 


Table  A3  Input  Data 


Description 

Symbol 

Vtilues 

first  principal  curvature  of  tool 

Kf 

0.0  (1/mm) 

second  principal  curvature  of  tool 

Kh 

0.0  (1/mm) 

first  principal  direction  of  tool 

e/ 

[0.9397  0.2962  0.1710]^ 

second  principal  direction  of  tool 

ek 

[0.0  -  0.5  0.8660]^ 

angular  velocity  of  tool  (1/sec) 

uW 

[0.0  0.0  0.0]^ 

angular  velocity  of  gear  (1/sec) 

[0.0  0.0  -0.2003f 

transfer  velocity  of  tool 

(mm/sec) 

vi;> 

[0.0  58.6588  0.0]^ 

transfer  velocity  of  gear 

(mm/sec) 

vg> 

[-0.2096  58.7218  0.0]^ 

surface  normal  of  tangent  point 

n 

[-0.3420  0.8138  0.4698]^ 

Table  A4  Output  Data 


Description  (for  gear) 

Symbol 

Values 

first  principal  curvature 

-0.007836  (1/mm) 

second  principal  curvature 

0.0  (1/mm) 

first  principal  direction 

[0.9219  0.3874  0.0]^ 

second  principal  direction 

[0.1820  -  0.4331  0.8827]^ 
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Fig.  A1  Principal  Directions 
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Appendix  B.  Contact  Ellipse 


Determination  of  Dimensions  and  Orientation  of  Instantaneous  Contact  Ellipse 

The  gear  tooth  surfaces  are  in  point  contact  at  every  instant.  Due  to  elastic  deformation 
of  gear  tooth  surfaces  the  contact  is  spread  over  an  elliptical  area  and  the  center  of  the  ellipse 
coincides  with  the  instantaneous  contact  point.  The  bearing  contact  is  formed  as  the  set  of 
instantaneous  contact  ellipses. 

The  dimensions  and  orientation  of  the  instantaneous  contact  ellipse  can  be  determined 
iisiTig  the  data  about  the  principal  curvatures  and  directions  of  the  contacting  surfaces,  and 
the  elastic  approach  of  the  surfaces.  The  elastic  approach  depends  on  the  applied  load  but 
we  will  consider  it  as  a  given  value  that  is  known  from  experimental  data. 

The  determination  of  the  instantaneous  contact  ellipse  is  based  on  the  following  equations 
(proposed  by  Litvin): 


cos 


gi  —  g2  cos  2<t _ 

(^1  -  23^1^2  cos  2a  + 


m 


sin2a^^^  = 


g2  sin  2a 

{gl  -  2pi52  cos  2a  +  giyf^ 


(B2) 


s 

1/2 

s 

2a  =  2 

A 

,  26  =  2 

B 

1/2 


where 


(B3) 


A  =  -  /eg)  -  (gf  -  2^i^2  cos  2a  +  giy^^]  (B4) 

B  =  -  4^  +  {gl  -  2^152  cos  2a  +  gj^^]  {B5) 


_  .(*)  4.  ^(0 

/Cj.  —  Kj  -r  Kjj  , 


n-  -  -  k(’) 

g^  —  Kj  —  Kjj 


(B6) 
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Here  (fig.  Bl)  is  the  angle  that  is  formed  by  axis  rj  of  the  contact  ellipse  with  the  unit 
vector  ej  ^  of  the  principal  direction  on  surface  Si;  (T  is  the  angle  formed  by  unit  vectors 
Cj  ^  and  of  the  principal  directions  of  the  contacting  surfaces;  2a  and  26  are  the  axes  of 
the  contact  ellipse;  S  is  the  elastic  approach;  and  Kj  ^  and  Kjj  are  two  principal  directions  of 
tooth  surface  i. 

Numerical  Example 

The  input  and  output  for  the  determination  of  the  contact  ellipse  are  shown  in  Tables 
Bl  and  B2. 


Table  Bl  Input  Data 


Description 

Symbol 

Values 

pinion  first  principal  curvature 

-0.001543  (1/mm) 

pinion  second  principal  curvature 

0.03907  (l/mm) 

pinion  first  principal  direction 

e. 

[0.1751  -  0.4361  0.8827]^ 

pinion  second  principal  direction 

e? 

[0.9222  0.3865  0.0080]^ 

gear  first  principal  curvature 

Kf 

-0.007836  (l/mm) 

gear  second  principal  curvature 

Kh 

0.0  (l/mm) 

gear  first  principal  direction 

e/ 

[-0.9219  -  0.3874  O.Op 

gear  second  principal  direction 

[-0.1820  0.4331  0.8827]^ 

elastic  approach 

6 

0.007  (mm) 

Table  B2  Output  Data 


Description 

Symbol 

Values 

long  axis  of  contact  ellipse 

2a 

6.026  (mm) 

short  axis  of  contact  elhpse 

26 

1.092  (mm) 

angle  between  long  axis  and 

the  first  principal  direction  of  pinion 

ad) 

89.87  (deg) 
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Fig.  B1  Orientation  and  dimensions  of  contact  ellipse 
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Appendix  C.  Directions  for  Users  of  Application  of  Computer  Program 
C.l  Introduction 

The  name  of  the  computer  program  is  HELTCA.FOR.  It  is  written  in  FORTRAN?? 
language.  The  operating  system  is  CMS-9.0.  A  subroutine  DNEQNF  to  solve  a  system  of 
nonlinear  equations  should  be  available  in  the  Math-Library  or  working  environment.  The 
subroutine  is  not  included  in  the  program.  The  program  will  call  the  subroutine  DNEQNF 
several  times. 

C.2  Input  Block 

The  input  block  consists  of  three  parts:  (1)  design  parameters  of  pinion  and  gear;  (2)  the 
controlled  modification  parameters;  and  (3)  parameters  for  TCA. 

Part  i.  D&ign  parameters  of  pinion  and  gear 

In  the  beginning  of  the  computer  program,  you  can  read  the  following  lines: 

C  All . COEFFICIENT  FOR  TRANSFORMATION  OF  DEGREE  TO  RADIAN 

All=DACOS(-l.D0)/180.D0 

C 

C...  KHD=1  FOR  RIGHT-HAND  PINION  AND  LEFT-HAND  GEAR 
C...  KHD=2  FOR  LEFT-HAND  PINION  AND  RIGHT-HAND  GEAR 

KHD=2 

If  you  write  “KHD=1’’ ,  the  computer  will  use  the  necessary  equations  for  the  case  of 
right-hand  pinion  and  left-hand  gear.  The  computations  will  be  accomplished  for  a  left- 
hand  pinion  and  right-hand  gear  if  you  use  “KHD=2”. 

Then,  the  variable  definition  for  the  pinion  follows: 

C 

C...  INPUT  THE  DESIGN  PARAMETERS  OF  PINION 
C... 

C...  TNI . GEAR  NUMBER  OF  TEETH 

C...  PNl . NORMAL  DIAMETRAL  PITCH  (1/MM) 

C...  PSINl . NORMAL  PRESSURE  ANGLE  (RAD.) 

C...  BETAPl . HELICAL  ANGLE  ON  PITCH  CYLINDER  (RAD.) 

C...  ADGl . ADDENDUM  (MM) 


53 


c...  DEGl . DEDENDUM  (MM) 

C...  LAMDPl . LEAD  ANGLE  ON  PITCH  CYLINDER  (RAD.) 

C...  FWl . FACE  WIDTH  (MM) 

C...  RPTl . RADIUS  OF  PITCH  CYLINDER  (MM) 

C...  RBTl . RADIUS  OF  BASE  CYLINDER  (MM) 

C...  RATI . RADIUS  OF  ADDENDUM  CYLINDER  (MM) 

C...  RDTl . RADIUS  OF  DEDENDUM  CYLINDER  (MM) 

C...  PSITl . TRANSVERSE  PRESSURE  ANGLE  (RAD.) 

C...  LAMDB1....LEAD  ANGLE  ON  BASE  CYLINDER  (RAD.) 

C 

In  accordance  with  our  numerical  example  (see  Table  1),  the  following  data  would  be 


used: 


TN1=20.D0 

PN1=5.D0/25.4D0 

PSIN1=A11*20.D0 

LAMDP1=A11*60.D0 

BETAP1=A11*30.D0 

FW1=25.4D0*1.6D0 

ADG1=1.D0/PN1 

DEG1=1.25D0/PN1 


The  following  variables  are  used  for  the  gear: 

C 

C...  INPUT  THE  DESIGN  PARAMETERS  OF  GEAR 
C 

C...  TN2 . GEAR  NUMBER  OF  TEETH 

C...  PN2 . NORMAL  DIAMETRAL  PITCH  (1/MM) 

C...  PSIN2 . NORMAL  PRESSURE  ANGLE  (RAD.) 

C...  BETAP2 . HELICAL  ANGLE  ON  PITCH  CYLINDER  (RAD.) 

C...  ADG2 . ADDENDUM  (MM) 

C...  DEG2 . DEDENDUM  (MM) 

C...  LAMDP2 . LEAD  ANGLE  ON  PITCH  CYLINDER  (RAD.) 

C...  FW2 . FACE  WIDTH  (MM) 

C  RPT2 . RADIUS  OF  PITCH  CYLINDER  (MM) 

C  ..  RBT2 . RADIUS  OF  BASE  CYLINDER  (MM) 

C  RAT2 . RADIUS  OF  ADDENDUM  CYLINDER  (MM) 

C  RDT2 . RADIUS  OF  DEDENDUM  CYLINDER  (MM) 

C...  PSIT2 . TRANSVERSE  PRESSURE  ANGLE  (RAD.) 

C...  LAMDB2 . LEAD  ANGLE  ON  BASE  CYLINDER  (RAD.) 

C 


From  the  design  parameters  listed  in  Table  1,  we  have 
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TN2=100.D0 

PN2=5.D0/25.4D0 

All=DACOS(-l.D0)/180.D0 

PSIN2=A11*20.D0 

LAMDP2=A11*60.D0 

BETAP2=A11*30.D0 

FW2=25.4D0*1.6D0 

ADG2=1.D0/PN2 

DEG2=1.25D0/PN2 


If  KHD=2,  the  computer  program  will  change  the  values  of  some  design  parameters  as 
follows 


IF(KHD.EQ.2)  THEN 
LAMDP1=-LAMDP1 
BETAP1=-BETAP1 
LAMDP2=-LAMDP2 
BETAP2=-BETAP2 
ENDIF 

The  computer  program  call  the  following  subroutines  “DATATl”  and  “DATAT2’’  to 
calculate  other  tooth  element  proportions  and  output  the  whole  data  in  file  55  (see  below). 

CALL  DATATl 
CALL  DATAT2 


Part  2.  Control  of  modification  parameter  for  application  in  the  TCA  program 
At  this  stage  we  can  read 


C  ...  THE  FOLLOWING  DATA  IS  FOR  THE  TO-BE  CONTROLLED  MODIFICATION 
PARAMETERS 


AA...  MODIFICATION  PARAMETER  OF  GEAR 

AP...  MODIFICATION  PARAMETER  OF  PINION  RACK-CUTTER 

DP  .  TANGENT  POINT  N  OF  PROFILES  OF  PINION  k  GEAR  RACK-CUTTERS 

THET2P..INITIAL  ANGLE  FOR  MODIFICATION  OF  GEAR  (RAD.) 


You  must  input  the  four  parameters  for  modification  of  pinion  and  gear  surfaces,  for 


example: 


AA=-0.0014D0 

THET2P=-0.08D0 

AP=-0.0008 

DP=-DSIN(PSIN1)*DCOS(-1.DO)*RPT1/TN1/8.DO 
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The  above  four  controlling  parameters  should  be  tried  several  times  is  order  to  obtain 
better  contact  pattern  and  transmission  errors  optimal  for  a  given  design. 


Part  3.  Parameters  for  TCA 

In  this  part,  the  alignment  errors  expected  should  be  input: 


C 

C  PARAMETERS  FOR  TCA 
C 

C..  KM  ..  SWITCH  1  FOR  CROSSING  ANGLE  A7,  &  2  FOR  INTERSECTION 
ANGLE  A7„ 

C..  DGAM..ANGLE  OF  MISALIGNMENT(CROSSING  OR  INTERSECTION)  (ARC 
MINUTE) 

C..  DEE...  CHANGE  OF  CENTER  DISTANCE  (MM) 

C 

If  a  crossing  angle  of  misalignment  is  considered,  should  input  “KM=1”.  Input  “KM=2” 
if  an  intersection  angle  of  misalignment  is  considered.  For  an  aligned  gear  drive,  input 
“DGAM=0.0’’ .  For  instance,  if  =  4',  AE  =  0,  input  the  following  lines: 


KM=1 
C  KM=2 

DGAM=4.0*A11/60.D0 

DEE=O.OOODO 

Then  you  will  read  the  following  sentence: 

C 

C....  THE  INPUT  BLOCK  IS  READY 

f 

Usually,  you  cannot  make  changes  anything  after  this  step. 


C.3  Output  Block 

After  the  input  block  is  filled  out,  you  can  read  the  following  explanation  for  the  output 


files: 

C 

C...  OUTPUT  DATA  FILES  ARE  THE  FOLLOWINGS 
C 

C...  FILE  55...  TOOTH  PROPORTIONS  OF  PINION  AND  GEAR 
C...  FILE  85...  CONTACT  PATH  ON  PINION  SURFACE  (2D) 

C...  FILE  86...  CONTACT  PATH  ON  GEAR  SURFACE  (2D) 

C...  FILE  87...  LENGTH  AND  DIRECTIONS  OF  CONTACT  ELLIPSE  ON  PINION 
AND  GEAR  SURFACES  (2D) 
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C...  FILE  90...  TRANSMISSION  ERRORS 
C 

C . 

C 

1)  File  55 

In  File  55  the  information  about  the  pinion  and  gear  are  listed. 

2)  Files  85  and  86 

There  are  two  coordinates  in  File  85  for  eadi  contact  point  of  the  pinion:  Radial  (Xp+yp)*^’* 
and  axial  Zp  .  There  are  two  coordinates  in  File  86  for  each  contact  point  of  the  gear:  Radial 

(^3  +  2/3)°'® 

3)  File  87 

There  are  5  values  in  File  87  for  each  pair  of  contact  points  of  the  pinion  and  gear.  The 
first  one  is  the  value  of  the  major  semi-axis  of  the  contact  ellipse.  The  second  and  third 
values  are  the  cosine  directions  of  the  major  axis  of  the  contact  ellipse  on  the  pinion  tooth 
surface.  The  last  two  values  are  the  cosine  direction  of  the  major  axis  of  the  contact  ellipse 
on  the  gear  tooth  surface. 

C.4  Computer  Program 
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c 


TCA  FOR  MODIFIED  HELICAL  GEARS . 

PROGRAM  HELTCA 
IMPLICIT  REAL*8(A“H,0-Z) 

REAL*8  XI(9) ,X(9) ,F(9) 

REAL*8  LF1(3,3) ,LH2(3,3) ,L1F{3,3) ,L2H(3,3) 

REAL*8  R1F(3) ,R2F(3) ,N1F(3) ,N2F(3) 

R£AL*8  DPHI(2,180) ,DDPHI(2,180,4) 

REAL*8  ELAL(180) , ELI (2, 180) ,EL2 (2,180) 

REAL*8  LFH(3,3) ,LHF(3,3) 

REAL*8  Vl(3) ,V2(3) ,V3(3) ,V4(3) ,V5(3) 

REAL*8  RG2(3) ,NG2(3) 

REAL*8  RG1(3) ,NG1(3) 

REAL*8  UI(3) ,DJ(3) ,UK(3) 

REAL*8  LAMDP2,LAHDB2 
REAL*8  LAHDP1,LAMDB1 
REAL*8  DT2,KT21,KT22 

REAL*8  AVC1(3) ,VTR1(3) ,AVC2(3) ,VTR2(3) 

REAL*8  KSIG1,KSI62,KFF,KHH 

REAL*8  EFN(3)  ,EHN(3)  ,W1VT2(3)  ,WV12(3)  ,W2VT1(3)  ,KF,KH,KS,KQ 
R£AL*8  KM2,KT2 
REAL*8  KFP,KHP,KFG,KHG 
COMMON  /A300/  ES(3),EQ(3) 

COMMON  /A310/  KFF,KHH 
COMMON  /A340/  EFF(3) ,EHH(3) 

COMMON  /A360/  A,B,SI(3) ,FI(3) 

COMMON/A200/  Wl(3) ,W2(3) ,W12(3) ,VT1(3) ,VT2(3) ,V12(3) 

COMMON  /A210/  EX(3) ,EF(3) ,EH(3) 

COMMON  /A212/  EF2 (3) ,EH2 (3) ,KF2 ,KH2 
COMMON  /A220/  KF,KH 
COMMON  /A230/  ET(3),EM(3) 

COMMON  /A380/  KS,KQ 

COMMON  /A400/  VT11(3) ,VT12 (3) ,VT21(3) ,VT22 (3) 

COMMON  /A401/  KHP , KFP , KHG , KFG 
EXTERNAL  FCNG,FCNC,FCNT 
COMMON  /AXIS/  UI,UJ,UK 
COMMON  /NET/  RR,DD 
COMMON  /DD/  DF,KPRI 

COMMON  / DATT2 /  TN2 , PN2 , PS IN2 , BETAP2 , ADG2 , DEG2 , LAMDP2 , 

&  UP2 , FW2 , RPT2 , RBT2 , RAT2 , RDT2 , PSIT2 , LAMDB2 
COMMON  /DATT/  TNI , PNl , PSINl , BETAPl , ADGl , DEGl , LAMDPl , 

&  FW1,RPT1,RBT1,RAT1,RDT1,PSIT1,LAMDB1 
COMMON  /B2/XMP1 , YNPl , ZNPl , XNP2 , YNP2 , ZNP2 
COMMON  /B4/X1,Y1,Z1,XN1,YN1,ZN1 
COMMON  / B5 /X2 , Y2 , Z2 , XN2 , YN2 , ZN2 
COMMON  /B6/THET2P,DGPHI2 

COMMON  /Wl/  ETAW1,UPP1,SPP1,ETAW2,UPP2,SPP2 
COMMON  /SGI/  RG1,NG1,AP,DP 
COMMON  /SG2/  RG2,NG2,AA 

COMMON  /ATT/  PHI1,PHI2,R1F,R2F,N1F,N2F,LFH,CC,DGAM,DPHI2 
COMMON  /MVT/  LF1,LH2 
COMMON  /AST/  ICONT 
COMMON  /ATS/  DPHIl 

C...  All . COEFFICIENT  FOR  TRANSFORMATION  OF  DEGREE  TO  RADIAN 

All=DACOS(-l.D0)/180.D0 

C...  KHD=1  FOR  RIGHT-HAND  PINION  AND  LEFT-HAND  GEAR 
C...  KHD=2  FOR  LEFT-HAND  PINION  AND  RIGHT-HAND  GEAR 
KHD=2 

C...  INPUT  THE  DESIGN  PARAMETERS  OF  PINION 


C...  TNI . GEAR  NUMBER  OF  TEETH 

C...  PNl . NORMAL  DIAMETRAL  PITCH  (MM) 
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C...  PSINl . NORMAL  PRESSURE  ANGLE  (RAD.) 

C...  BETAPl . LEADING  ANGLE  OF  THE  HELIX  ON  PITCH  CYLINDER  (RAD.) 

C...  ADGl . ADDENDUM  (MM) 

C...  DEGl . DEDENDUM  (MM) 

C...  LAMDPl _ HELIX  ANGLE  ON  PITCH  CYLINDER  (RAD.) 

C...  FWl . FACE  WIDTH  (MM) 

C...  RPTl . RADIUS  OF  PITCH  CYLINDER  (MM) 

C...  RBTl . RADIUS  OF  BASE  CYLINDER  (MM) 

C...  RATI . RADIUS  OF  ADDENDUM  CYLINDER  (MM) 

C...  RDTl . RADIUS  OF  DEDENDUM  CYLINDER  (MM) 

C...  PSITl . TRANSVERSE  PRESSURE  ANGLE  (RAD.) 


C...  LAMDBl _ HELIX  ANGLE  ON  BASE  CYLINDER  (RAD.) 

TN1*20.D0 
PN1=5. DO/25. 4D0 
PSIN1=A11*20.D0 
LAMDP1=A11*60 . DO 
BETAP1=A11*3  0 . DO 
FW1=25.4D0*1.6D0 
ADG1=1.D0/PN1 
DEG1=1.25D0/PN1 

C...  INPUT  THE  DESIGN  PARAMETERS  OF  GEAR 


C. . .  TN2 . GEAR  NUMBER  OF  TEETH 

C. ..  PN2 . NORMAL  DIAMETRAL  PITCH  (1/MM) 

C. ..  PSIN2 . NORMAL  PRESSURE  ANGLE  (RAD.) 

C...  BETAP2 . LEADING  ANGLE  OF  THE  HELIX  ON  PITCH  CYLINDER  (RAD.) 

C. ..  ADG2 . ADDENDUM  (MM) 

C...  DEG2 . DEDENDUM  (MM) 

C. ..  LAMDP2 . HELIX  ANGLE  ON  PITCH  CYLINDER  (RAD.) 

C...  FW2i . FACE  WIDTH  (MM) 

C...  RPT2 . RADIUS  OF  PITCH  CYLINDER  (MM) 

C. ..  RBT2 . RADIUS  OF  BASE  CYLINDER  (MM) 

C. ..  RAT2 . RADIUS  OF  ADDENDUM  CYLINDER  (MM) 

C...  RDT2 . RADIUS  OF  DEDENDUM  CYLINDER  (MM) 

C. ..  PSIT2 . TRANSVERSE  PRESSURE  ANGLE  (RAD.) 

C...  LAMDB2 . HELIX  ANGLE  ON  BASE  CYLINDER  (RAD.) 


TN2*=100.D0 
PN2=5.D0/25.4D0 
All=DACOS(-l.D0)/180.D0 
PSIN2=A11*20.D0 
LAMDP2=A11*60 . DO 
BETAP2=A11* 3 0 . DO 
FW2=25.4D0*1.6D0 
ADG2=1.D0/PN2 
DEG2=1.25D0/PN2 
IF(KHD.EQ.2)  THEN 
LAMDP1=-LAMDP1 
BETAP1=-BETAP1 
LAMDP2=-LAMDP2 
BETAP2=-BETAP2 
ENDIF 

CALL  DATATl 
CALL  DATAT2 

C. . .  THE  FOLLOWING  IS  FOR  CONTROLLING  PARAMETERS 

C. . .  AA. . .  MODIFICATION  PARAMETER  OF  GEAR 

C.  ..  AP...  MODIFICATION  PARAMETER  OF  PINION  RACK-CUTTER 

C. .  DP.  .  .  .TANGENT  POINT  N  OF  PROFILES  OF  PINION  &  GEAR  RACK-CUTTERS 

C..  THET2P. .INITIAL  ANGLE  FOR  MODIFICATION  OF  GEAR  (RAD.) 

C 

C  PARAMETERS  FOR  TCA 
C 
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C. .  KM  . .  SWITCH  1  FOR  CROSSING  &  2  FOR  INTERSECTING  MISALIGNMENT 
C..  DGAM.  .ANGLE  OF  DISALIGNMENT (CSOSSING  OR  INTERSECTING)  (ARC  MINUTE) 
C..  DEE...  CHANGE  OF  CENTER  DISTANCE  (MM) 

KM=1 
C  KM=2 

DGAH^O.DO 
DEE-0. OOODO 
AA=-0.0014DO 
THET2P*-0.08D0 
AP=-0.0008 

DP“-DSIN(PSIN1) *DCOS(-1.DO) *RPTl/TNl/8 .  DO 
C...  THE  INPUT  BLOCK  IS  OVER  HERE 
C 

C...  OUPTPUT  DATA  FILES  ARE  THE  FOLLOWINGS 
C 

C...  FILE  55...  TOOTH  PROPORTIONS  OF  PINION  AND  GEAR 
C...  FILE  85...  CONTACT  PATH  ON  PINION  SURFACE  (2D) 

C...  FILE  86...  CONTACT  PATH  ON  GEAR  SURFACE  (2D) 

C...  FILE  87...  DIRECTIONS  OF  LONG  AXIS  OF  CONTACT  ELLIPS  (2D) 

C...  FILE  90  ...TRANSMISSION  ERRORS 
DO  901  1=1,3 
UI(I)=0.D0 
UJ(I)=0.D0 
DK(I)=0.D0 
901  CONTINUE 
UI(1)=1.D0 
UJ(2)=1.D0 
UK(3)=1.D0 

C  ..  EE2.....  GEARRATIO 
EE2=TN1/TN2 

C  . .  CC .  CENTER  DISTANCE  OF  GEAR  DRIVE 

CC=RPT1+RPT2+DEE 

C  . .  CALCULATE  CONTACT  POINT  ON  MEAN  SECTION  WITHOUT  MISALIGNMENT 
ICONT=l 
N=6 

ERRREL=0 . lD-6 
ITMAX=1000 
C 

CALL  MIAL(KM,0.D0,LFH) 

C 

XI(1)=0.D0 

XI(2)=O.DO 

XI(3)=0.D0 

XI(4)=0.D0 

XI(5)=O.DO 

XI(6)=0.D0 

DD=0.D0 

C 

CALL  DNEQNF(FCNG,ERRREL,N,ITMAX,XI,X,FNORM) 

C 

PHISS1=X(6) 

PHI1MEA=X(6) 

C  ..  CALCULATE  CONTACT  POINT  ON  EDGE  SECTION  WITHOUT  MISALIGNMENT 
DD=-0 . 5D0*FW1 
ERRREL=0 . lD-8 
c 

CALL  DNEQNF(FCNG,ERRREL,N,ITMAX,XI,X,FN0RM) 

c 

PHISS2=X ( 6 ) 

C  ..  CALCULATE  CONTACT  POINT  ON  EDGE  SECTION  WITH  MISALIGNMENT 
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N=6 

ERRREIr=0.1D-6 

ITMAX=1000 

C 

CALL  HIAL(KH,D6AM,LFH) 

C 

XI(1)=X(1) 

XI(2)=X(2) 

XI(3)=X(3) 

XI(4)=X(4) 

XI(5)=X(5) 

XI(6)=X(6) 

C 

CALL  DNEQNF{FCNG,ERRREL,N,ITMAX,XI,X,FNORM) 

C 

1330  CONTINUE 
PHI1=X{6) 

PHI1STA«X(6) 

C  . .  THE  FOLLOWING  IS  FOR  TCA 
ICONT=2 
N=5 

EERREL=0.1D-5 

ITMAX=400 

STP=(PHISS2-PHISS1)/36.D0 

NN=72 

DO  1010  1=1, NN 

PHI1=PHI1STA-(I-1) *STP 
XI(1)=X(1) 

XI(2)=X(2) 

XI(3)=X(3) 

XI(4)=X(4) 

XI(5)=X(5) 

C 

CALL  DNEQNF(FCNG,ERRREL,N,ITMAX,XI,X,FNORM) 

IF(RG1(3) .GT. (0.5D0*FW1))  GOTO  1011 

C.. .  ROTATING  VELOCITY  OF  CUTTER  IN  CUTTER  SYSTEM 
W1(1)=0.D0 
W1(2)=0.D0 
W1(3)=0.D0 

C...  ROTATING  VELOCITY  OF  PINION  IN  CUTTER  SYSTEM 
W2(1)=0.D0 
W2(2)=0.D0 
W2(3)=1.D0 

C...  NORMAL  OF  CUUTER  IN  CUTTER  SYSTEM 
EX{1)=XNP1 
EX(2)=YNP1 
EX(3)=ZNP1 

C. . .  TRANSFER  VELOCITIES  OF  CUTTER  AND  PINION  IN  CUTTER  SYSTEM 
CALL  EQVEC(VT1,VT11) 

CALL  EQVEC(VT2,VT12) 

C. . .  RELATIVE  VELOCITIS  OF  CUTTER  WRT.  PINION  IN  CUTTER  SYSTEM 
CALL  ADDVEC(W12,W1,W2,--1.D0) 

CALL  ADDVEC(V12,VT1,VT2,-1.D0) 

KF=KFP 

KH=KHP 

C...  PRINCIPAL  CURVATUTES  AND  DIRECTIONS  OF  PINION 
CALL  CURVT(l) 

C...  PRINCIPAL  DIRECTIONS  OF  PINION  IN  PINION  SYSTEM 
CALL  EQVEC(EFF,ES) 
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CALL  EQVEC(EHH,EQ) 

C. . .  ROTATING  VELOCITY  OF  CUTTER  IN  CUTTER  SYSTEM 
W1(1)-0.D0 
W1{2)=0.D0 
W1(3)=0.D0 

C...  ROTATING  VELOCITY  OF  GEAR  IN  CUTTER  SYSTEM 
W2(1)=0.D0 
W2(2)=0.D0 
W2(3)=“DGPHI2 

C. . .  NORMAL  OF  CUUTER  IN  CUTTER  SYSTEM 
EX(1)=XNP2 
EX(2)=YNP2 
EX(3)=ZNP2 

C. . .  TRANSFER  VELOCITIES  OF  CUTTER  AND  PINION 
CALL  EQVEC(VT1,VT21) 

CALL  EQVEC(VT2,VT22) 

CALL  EQVEC(EF,EF2) 

CALL  EQVEC(EH,EH2) 

KF=KF2 

KH=KH2 

C...  RELATIVE  VELOCITIS  OF  CUTTER  WRT.  GEAR  IN  CUTTER  SYSTEM 
CALL  ADDVEC(W12,W1,W2,-1.D0) 

CALL  ADDVEC(V12,VT1,VT2,-1.D0) 

KF=KFG 

KH=KHG 

C. ..  PRINCIPAL  CURVATUTES  AND  DIRECTIONS  OF  PINION 
CALL  CURVT(2) 

C...  PRINCIPAL  DIRECTIONS  OF  PINION  IN  PINION  SYSTEM 
CALL  TRANSM(L1F,LF1) 

CALL  MAVEC(V1,LH2,ES) 

CALL  MAVEC(V2,LFH,V1) 

CALL  MAVEC(ES,L1F,V2) 

CAT.L  MAVEC(V1,LH2,EQ) 

CALL  MAVEC(V2,LFH,V1) 

CALL  MAVEC(EQ,L1F,V2) 

C...  CONTACT  ELLIPSES 
CALL  ELLIP 
ELAL(I)=A 

C...  AXIS  OF  PINION  ELLIPSE  ON  TANGENT  PLANE 
CMSL  EQVEC(  VI, SI) 

CALL  DOTVEC(VlN,Vl,Vl) 

EL1(1,I)=DSQRT(V1(1)**2+V1(2)**2/V1N) 

ELI ( 2 , I ) =V1 ( 3 ) /DSQRT ( VIN) 

C...  AXIS  OF  GEAR  ELLIPSE  ON  TANGENT  PLANE 
CALL  TRANSM(LHF,LFH) 

CALL  TRANSM(L2H,LH2) 

CALL  MAVEC(V2,LF1,V1) 

CALL  MAVEC(V1,LHF,V2) 
rALL  MAVEC(V2,L2H,V1) 

CALL  EQVEC(V1,V2) 
rAT.L  D0TVEC(V1N,V1,V1) 

EL2 ( 1 , I) =-DSQRT (VI ( 1) **2+Vl ( 2 ) **2 /VIN) 

EL2 ( 2 , I ) =V1 ( 3 ) /DSQRT ( VIN) 

AGl^DATAN (ELI (2 , 1) /ELI (1 , 1) )  /All 

AG2=DATAN { EL2 ( 2 , I ) / EL2 ( 1 , I ) ) / All 

DPHI(1,I)»PHI1 

DPHI (2 , I) =PHI2-TN1/TN2*PHI1 

DDR1=DSQRT (RGl (1) **2+RGl (2 ) **2 ) 

DDR2=DSQRT (RG2 (1) **2+RG2 (2 ) **2 ) -210 . DO 
KM=4 
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non 


0U1=FI.0AT  (I-l)  /  FLOAT  (KM) 

0UP=AINT(0U1) 

IF(OUl.EQ.OUP)  THEN 

WRITE (85,*)  DDR1,RG1(3) 

WRITE (86,*)  DDR2,RG2(3) 

C  WRITE(87,*)  EliAL(I)  ,EL1(1,I)  ,EL1(2,I)  ,EL2(1,I)  ,EL2(2,I) 

ENDIF 

1010  CONTINUE 

1011  CONTINUE 
AD=-10000.DO 

DO  1090  I=1,NN,KM 

IF(DPHI(2,I) .NE.O.DO)  THEN 
IF(DPHI(2,I) .GT.AD)  THEN 
AD=DPHI(2,I) 

ENDIF 

ENDIF 

1090  CONTINUE 

SS=360.D0/TN1 
DO  1020  J=l,4 

DO  1030  I=1,NN,KM 

BD=(DPHI(2,I)-AD)/A11*3600.DO 
DDPHI ( 1 , I , J) “DPHI (1,1) / All+SS* ( J-1 ) 

DDPHI(2,I,J)=BD 
1030  CONTINUE 
1020  CONTINUE 
KM=4 

DO  1040  J=l,4 

DO  1050  I=1,NN,KM 

IF(DDPHI(1, I, J) .NE.O.DO)  THEN 

WRITE(90,*)  DDPHI(1,I,J) ,DDPHI(2,I, J) 

ENDIF 

1050  CONTINUE 
1040  CONTINUE 
KM=4 

DO  1060  I=1,NN,KM 

WRITE(87,*)  ELAL(I) ,EL1{1,I)  , ELI (2 , 1) , EL2 (1 , 1) ,  EL2  (2 , 1) 
1060  CONTINUE 

WRITE (6,*)  ******  PROGRAM  FINISHED  ******* 

STOP 
END 

...THE  SUBROUTINE  IS  FOR  TCA 

SUBROUTINE  FCNG(X,F,N) 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL* 8  X(N),F(N) 

REAL*8  LF1(3,3) ,LH2(3,3) ,LFH(3,3) 

REAL*8  R1F(3) ,R2F(3) ,N1F(3) ,N2F(3) ,CO(3) 

REAL*8  AVC1(3) ,VTR1(3) ,AVC2(3)  ,VTR2(3) 

REAL*8  UI(3) ,UJ(3) ,UK(3) ,V1(3)  ,V2(3)  ,V3(3) ,V4(3) ,V5(3) 
REAL*8  UP,LAMDP,LAMDB 
REAL*8  KF2,KH2 
REAL*8  KFP,KHP,KFG,KHG 
COMMON  /AXIS/  UI,UJ,UK 
COMMON  /SGI/  RG1,NG1,AP,DP 
COMMON  /SG2/  RG2,NG2,AA 

COMMON  /ATT/  PHI1,PHI2,R1F,R2F,N1F,N2F,LFH,CC,DGAM,DPHI2 
COMMON  /AST/  ICONT 
COMMON  /MVT/  LF1,LH2 
REAL*8  KSIG1,KSIG2,KFF,KHH 
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REAL*8  EFN(3)  ,EHN(3)  ,W1VT2{3)  ,WV12(3)  ,W2VT1(3)  ,KF,KH,KS,KQ 

REAL*8  KM2,KT2 

COMMON  /A300/  ES(3),EQ(3) 

COMMON  /A310/  KFF,KHH 
COMMON  /A340/  EFF(3) ,EHH(3) 

COMMON  /A340/  EFl (3) , EHl (3) 

COMMON  /A360/  A, B, SI (3) , FI(3) 

COMMON/A200/  Wl(3) ,W2(3) ,W12(3) ,VT1(3) ,VT2(3) ,V12(3) 

COMMON  /A210/  EX(3) ,EF(3) ,EH(3) 

COMMON  /A212/  EF2 (3) ,EH2 (3) ,KF2,KH2 
COMMON  /A220/  KF,KH 
COMMON  /A230/  ET(3),EM(3) 

COMMON  /A380/  KS,KQ 

COMMON  /A400/  VT11(3) ,VT12 (3) , VT21 (3) , VT22 (3) 

COMMON  /A401/  KHP,KFP,KHG,KFG 
REAL*8  RG1(3) ,N61(3) 

REAL*8  RG2(3) ,NG2(3) 

REAL*8  LAMDP2^LAMDB2 
REAL*8  LAMDPl.LAMDBl 

COMMON  /DATT2  /  TN2 ,  PN2  ,  PSIN2 ,  BETAP2  ,  ADG2 ,  DEG2  ,  LAMDP2 , 

&  UP2,FW2,RPT2,RBT2,RAT2,RDT2,PSIT2,LAMDB2 
COMMON  /B2/XNPl,YNPl,ZNPl,XNP2,yNP2,ZNP2 
COMMON  /B4/X1,Y1,Z1,XN1,YN1,2N1 
COMMON  /B5/X2,Y2,Z2,XN2,YN2,ZN2 
COMMON  /B6/THET2P,DGPHI2 

COMMON  /Wl/  ETAW1,UPP1,SPP1,ETAW2,UPP2,SPP2 
COMMON  /NET/  RR,DD 

COMMON  /DATT/  TN1,PN1,PSIN1,BETAP1,ADG1,DEG1,LAMDP1, 

&  FW1,RPT1,RBT1,RAT1,RDT1,PSIT1,LAMDB1 
UPP1=X{1) 

ETAW1»X(2) 

UPP2=X(3) 

ETAW2=X(4) 

PHI2=X{5) 

IF(ICONT.EQ.l)  THEN 
PHI1=X(6) 

ENDIF 

RPT=RPT1 

CNST=DARCOS(-1.0D00) / 18 0.0 
PI=DARCOS (-1. ODOO) 

PI2=2.D0*PI 
PSIT=  PSITl 
PSIN=  PSINl 
BETAP=  BETAPl 
DSIN1=DSIN (ETAWl ) 

DCOSl=DCOS ( ETAWl ) 

ANFl=DARCOS (RBTl/RPTl) 

CINV1=DTAN(ANF1) -ANFl 
ANG1=PI/2.D0/TN1 
AM=PI/PN1/4.D0 
SA=AM*DCOS (BETAP) 

SS1=  RPTl* (ETAWl) 

DSS1=  RPTl 

C...  EQUATION  OF  MESHING  OF  PINION 
FF1=SPP1*DSIN (BETAP) 

FF2*=(-{UPP1+DP)-2.D0*AP**2*UPP1**3) 

FF3=DSIN  (PSIN) +2 . DO*AP*UPP1*DCOS (PSIN) 

F5=FF2/FF3*DCOS (BETAP) -SA-RPT1*ETAW1 
SPP1=  F5/DSIN (BETAP) 

C...  SURFACE  OF  PINION  RACK  CUTTER 
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XP= (UPPl+DP) *DCOS (PSIN) 

YP= ( (UPPl+DP) *DSIN (PSIN) +AM) *DCOS (BETAP) + 

&  AP*UPPl**2*DCOS (PSIN) *DCOS (BETAP) +SPP1*DSIN (BETAP) 

ZP=- ( (UPPl+DP) *DSIN (PSIN) +AM) *DSIN (BETAP) - 
&  AP*UPP1**2*DC0S (PSIN) *DSIN (BETAP) +SPP1*DC0S (BETAP) 

Xl=  DCOS (ETAWl) *XP+DSIN (ETAWl) *YP+RPTl*DCOS (ETAWl) 

&+SSl*DSIN (ETAWl) 

Y1=-DSIN (ETAWl) *XP+DCOS (ETAWl) *YP-RPT1*DSIN (ETAWl) 

&+SSl*DCOS (ETAWl) 

Z1=ZP 

NORMAL  OF  RACK  CUTTER 
PNN=DSQRT(1.D0+(2.D0*AP*UPP1) **2) 

XNPP=-2 .D0*AP*UPP1/PNN 

YNPP=1.D0/PNN 

ZNPP®O.DO 

XNP1= (DCOS (PSIN) *XNPP-DSIN (PSIN) *YNPP) 

YNP1=  (DSIN (PSIN) *DCOS (BETAP) *XNPP+DCOS (PSIN) *DCOS (BETAP) *YNPP) 
ZNP1=(-DSIN(PSIN)  *DSIN(BETAP)  *XNPP-DCOS  (PSIN)  *DSIN(BETAP)  *YNPP) 
NORMAL  OF  PINION  IN  SI 
XN1=  DCOS (ETAWl) *XNP1+DSIN( ETAWl) *YNP1 
YN1=“DSIN (ETAWl) *XNPl+DCOS (ETAWl) *YNP1 
ZN1=  ZNPl 

DX1=-DSIN (ETAWl) *XP+DCOS (ETAWl) *YP+SSl*DCOS( ETAWl) 

DYl»-DCOS (ETAWl) *XP-DSIN (ETAWl) *YP-SS1*DSIN (ETAWl) 

DZ1=0.D0 

RG1(1)=X1 

RG1(2)=Y1 

RG1(3)*Z1 

NG1(1)=XN1 

NG1(2)=YN1 

NG1(3)»ZN1 

KHP=O.DO 

KFP=2 . D0*AP/ (DSQRT ( ( 1 . D0+ ( 2 . D0*AP*UPP1) **2 ) ) **3 ) 

EFl (1) =DCOS (PSIN) 

EFl ( 2 ) =DSIN (PSIN) *DCOS ( BETAP) 

EF1(3)=--DSIN(PSIN)  *DSIN  (BETAP) 

EH1(1)=0.D0 
EH1(2)=  DSIN (BETAP) 

EH1(3)=  DCOS (BETAP) 

EF ( 1) =  DCOS (ETAWl) *EH1 ( 1) +DSIN (ETAWl) *EH1 (2 ) 

EF ( 2 ) =-DSIN ( ETAWl) *EH1 ( 1 ) +DCOS ( ETAWl) *EH1 ( 2 ) 

EF(3)=  EH1(3) 

EH ( 1 ) =  DCOS ( ETAWl ) *EF1 ( 1 ) +DSIN ( ETAWl ) *EF1 ( 2 ) 

EH ( 2 ) =-DSIN ( ETAWl) *EF1 ( 1 ) +DCOS (ETAWl ) *EF1 ( 2 ) 

EH(3)=  EF1(3) 

VELOCITY  OF  RACK-CUTTER  IN  SP 
VT11(1)=0.D0 
VT11(2)*RPT1 
VT11(3)=0.D0 

VELOCITY  OF  PINION  IN  SP 
VT12(1)=-(YP+SS1) 

VT12(2)=  XP+RPTl 

VT12(3)=0.D0 

RPT=RPT2 

PSIT=  PSIT2 

PSIN=  PSIN2 

BETAP=  BETAP2 

SS2=  RPT2*(ETAW2/TN2*TN1) 

DSS2*  RPT2/TN2*TN1 
TTD=AA*  (ETAW2-THET2P) **2 
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DTP=2 .  DO*AA*  ( ETAW2  -THET2P ) 

GPHI2=  ETAW2/TN2*TN1-TTD 
DGPHI2=  l.D0/TN2*TNl-DTP 
C...  EQUATION  OF  MESHING  OF  GEAR 

GG1=SPP2*DSIN (PSIN2) *DSIN (BETAP2) 

GG1*DSIN (PSIN2) *DSIN (BETAP2 ) 

GG2=(UPP2+DP) *DCOS (BETAP2) 

GG3=  RPT2*TNi/TN2*ETAW2*DSIN(PSIN2) 

GG4=SA*DSIN (PSIN2 ) 

GG5=  2  .D0*AA*TN2*  (ETAW2-THET2P)  *RPT2*DCOS  (PSIN2)  *DCOS  (BETAP2 ) 
GG6=  TN1-2.D0*AA*TN2*(ETAW2-THET2P) 

SPP2=- (GG2+GG3+GG4+GG5/GG6) /GGl 
C. . .  SURFACE  OF  GEAR  RACK  CUTTER 
XP= (UPP2+DP) *DCOS (PSIN) 

yP=  (  (UPP2+DP)  *DSIN  (PSIN)  +AM)  *DCOS  (BETAP)  +SPP2  *DSIN  (BETAP) 
ZP=-(  (UPP2+DP)  *DSIN(PSIN)  +AM)  *DSIN (BETAP)  +SPP2*DCOS  (BETAP) 
X2=-DCOS  (GPHI2 )  *XP+DSIN  (GPHI2  )  *yP+RPT2  *DCOS  (GPHI2 ) 
&+SS2*DSIN(GPHI2) 

Y2=-DSIN(GPHI2)  *XP-DCOS  (GPHI2)  *yP+RPT2*DSIN(GPHI2) 

&-SS2*DCOS (GPHI2) 

Z2=ZP 

C. . .  NORMAL  OF  RACK  COTTER  IN  SW 
XNP2=DSIN(PSIN) 
yNP2=-DCOS (PSIN) *DCOS (BETAP) 

ZNP2=  DCOS (PSIN) *DSIN (BETAP) 

XNP2=-DSIN(PSIN) 

yNP2=  DCOS (PSIN) *DCOS ( BETAP ) 

ZNP2=-DCOS (PSIN) *DSIN (BETAP) 

C  ..  NORMAL  OF  GEAR  IN  S2 

XN2=-DCOS (GPHI2) *XNP2+DSIN(GPHI2) *yNP2 
yN2=-DSIN(GPHI2) *XNP2-DCOS (GPHI2) *yNP2 
ZN2=  ZNP2 

DX2=(DSIN(GPHI2)  *XP+DCOS  (GPHI2)  *yp-RPT2*DSIN(GPHI2)  + 
&SS2*DCOS (GPHI2) ) *DGPHI2+DSS2*DSIN(GPHI2) 

Dy2=(-DCOS  (GPHI2)  *XP+DSIN(GPHI2)  *yP+RPT2*DCOS  (GPHI2)  + 
&SS2*DSIN(GPHI2) ) *DGPHI2-DSS2*DCOS (GPHI2) 

DZ2-0.D0 

RG2(1)=X2 

RG2(2)=y2 

RG2(3)=Z2 

NG2(1)=XN2 

NG2(2)=yN2 

NG2(3)=ZN2 

KFG=O.DO 

KHG=0 . DO 

EF2(l)=DCOS(PSIN) 

EF2 (2 ) =DSIN(PSIN) *DCOS (BETAP) 

EF2 (3 ) =-DSIN (PSIN) *DSIN (BETAP) 

EH2(1)=0.D0 
EH2(2)=  DSIN (BETAP) 

EH2(3)=  DCOS (BETAP) 

C  VELOCITY  OF  RACK-CUTTER  IN  SG 

VT21(1)=O.DO 
VT21(2)=®  RPTl 
VT21(3)=0.D0 
C  VELOCITY  OF  GEAR 

VT22 ( 1) =- ( YP+SS2 ) *DGPHI2 
VT22(2)=  (XP+RPT2)*DGPHI2 
VT22(3)=0.D0 
S1=DSIN(PHI1) 
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C1=DC0S(PHI1) 

S2»DSIN(PHI2) 

C2=DC0S(PHI2) 

LF1(1,1)=C1 

LF1(1,2)=S1 

LF1(2,1)=-S1 

LF1(2,2)=C1 

LF1(1,3)=0.D0 

LF1(2,3)=0.D0 

LF1(3,1)=0.D0 

LF1(3,2)*O.DO 

LF1(3,3)=1.D0 

LH2(1,1)=»C2 

LH2(1,2)=-'S2 

LH2(2,1)=S2 

LH2(2,2)=C2 

LH2(1,3)=0.D0 

LH2(2,3)=0.D0 

LH2(3,1)=0.D0 

LH2(3,2)=0.D0 

LH2(3,3)=1.D0 

CALL  MAVEC(V1,LH2,NG2) 

CALL  MAVEC(N2F,LFH,V1) 

CALL  MAVEC(V2,LH2,RG2) 

CALL  MAVEC(V3,LFH,V2) 

CALL  ADDVEC(R2F,V3,UI,CC) 

CALL  MAVEC(N1F,LF1,NG1) 

CALL  MAVEC(R1F,LF1,RG1) 

F(1)=(R1F{1)-R2F{1)) 

F(2)=(R1F(2)-R2F(2) ) 

F(3)=(R1F(3)-R2F(3)) 

F(4)=N1F(1)-N2F(1) 

F(5)*N1F(2)-N2F(2) 

IF(ICONT.EQ.l)  THEN 
F(6)=Z1-DD 
ENDIF 
RETURN 
END 

C. .  . 

c!!!  FOR  PINION  DATA  ** 

C.  •  • 

SUBROUTINE  DATATl 
IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  UP,LAMDP1,LAMDB1 

COMMON  /DATT/  TN1,PN1,PSIN1,BETAP1,ADG1,DEG1,LAMDP1, 

&  FW1,RPT1,RBT1,RAT1,RDT1,PSIT1,LAMDB1 
RPT1=TN1/ (2 . DO*PN1*DCOS (BETAPl) ) 

PSIT1=DATAN(DTAN(PSIN1) /DCOS (BETAPl) ) 

RBT1=RPT1*DC0S (PSITl) 

RAT1=RPT1+ADG1 

RDT1=RPT1-DEG1 

LAMDB1=DATAN(DTAN(LAMDP1) /DCOS (PSITl) ) 

BETAB1=DATAN (DTAN (BETAPl) *DCOS (PSITl) ) 

WRITE (5 5, 110) 

110  FORMAT(/2Xj  ^  ******************************************  ^  j 
h  2X,'*  DATA  OF  PINION  *',/ 

&  2X, ' ****************************************** , /) 

WRITE  (55 ,120)  TNI ,  PNl ,  PSINl ,  BETAPl ,  ADGl ,  DEGl ,  LAMDPl ,  FWl , 

&  RPT1,RBT1, RATI, RDTl, PSITl, lAMDBl 
120  FORMAT (2X, 'GEAR  NUMBER  OF  TEETH  TM1=',F14.7/ 
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&2X, 'NORMAL  DIAMETRAL  PITCH  (1/MM)  PN1=',F14.7/ 

&2X, 'NORMAL  PRESSURE  ANGLE  (RAD. )  PSIN1=' , F14 .7/ 

&2X, 'LEADING  ANGLE  OF  HELIX'  / 

&4X,'ON  PITCH  CYLINDER  (RAD.)  BETAP1= ' , F14 . 7 / 

&2X,  'ADDENDUM  (MM)  ADG1=' , F14 .7/ 

&2X, 'DEDENDUM  (MM)  DEG1=' , F14 .7/ 

&2X, 'HELIX  ANGLE  ON  PITCH  CYLINDER(RAD)  LAMDP1=' ,F14.7/ 
&2X, 'FACE  WIDTH  (MM)  FW1=',F14.7/ 

&2X, 'RADIUS  OF  PITCH  CYLINDER  (MM)  RPT1= ' , F14 . 7 / 

&2X, 'RADIUS  OF  BASE  CYLINDER  (MM)  RBT1=' ,  F14 . 7/ 

&2X, 'RADIUS  OF  ADDENDUM  CYLINDER  (MM)  RAT1= ' , F14 . 7 / 

&2X, 'RADIUS  OF  DEDENDUM  CYLINDER  (MM)  RDT1=' , F14 . 7/ 

&2X, 'TRANSVERSE  PRESSURE  ANGLE  (RAD>  PSIT1=' , F14 . 7/ 

&2X, 'HELIX  ANGLE  ON  BASE  CYLINDER  (RAD)  LAMDB1=' ,F14.7/) 

C.  . . 

RETURN 

END 

C..  . 

C  .  .  .  THE  SUBROUTINE  IS  FOR  DATA  OF  THEORITICAL  GEAR  SURFACE 

C. .  • 

SUBROUTINE  DATAT2 
IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  UP2,LAHDP2,LAMDB2 

COMMON  /DATT2/  TN2 ,PN2 ,PSIN2  ,BETAP2  ,  ADG2,DEG2 , LAMDP2 , 

&  UP2,FW2,RPT2,RBT2,RAT2,RDT2,PSIT2,LAMDB2 
RPT2=TN2/ (2.DO*PN2*DCOS(BETAP2) ) 

PSIT2=DATAN (DTAN (PSIN2) /DCOS (BETAP2 ) ) 

RBT2=RPT2*DCOS (PSIT2 ) 

RAT2=RPT2+ADG2 
RDT2=RPT2  -DEG2 

LAMDB2=DATAN(DTAN(LAMDP2) /DCOS(PSIT2) ) 

BETAB2=DATAN ( DTAN ( BETAP2 ) *DCOS ( PSIT2 ) ) 

DEL2*0.007 

C  WRITE(6,110) 

WRITE (55, 110) 

110  FORMAT (/2X, '*****************************************' , / 

&  2X,'*  DATA  OF  GEAR  2  *',/ 

&  2X, '****************************************** ^ /) 

WRITE  (55,120)  TN2 ,  PN2 ,  PSIN2 ,  BETAP2 , ADG2 ,  DEG2 ,  LAMDP2 , 

&  FW2 , RPT2 , RBT2 , RAT2 , RDT2 , PSIT2 , LAMDB2 , DEL2 

120  FORMAT (2X, 'GEAR  NUMBER  OF  TEETH  TN2=',F14.7/ 

&2X, 'NORMAL  DIAMETRAL  PITCH  (1/MM)  PN2=',F14.7/ 

&2X, 'NORMAL  PRESSURE  ANGLE  (RAD.)  PSIN2=' , F14 . 7/ 

&2X, 'LEADING  ANGLE  OF  HELIX'  / 

&4X,'ON  PITCH  CYLINDER  (RAD)  BETAP2=' , F14 . 7/ 

&2X, 'ADDENDUM  (MM)  ADG2=' , F14 . 7/ 

&2X, 'DEDENDUM  (MM)  DEG2=' , F14 . 7/ 

&2X, 'HELIX  ANGLE  ON  PITCH  CYLINDER(RAD)  LAMDP2=' , F14 . 7/ 
&2X, 'FACE  WIDTH  (MM)  FW2=',F14.7/ 

&2X, 'RADIUS  OF  PITCH  CYLINDER  (MM)  RPT2=' , F14 . 7/ 

&2X, 'RADIUS  OF  BASE  CYLINDER  (MM)  RBT2=' , F14 . 7/ 

&2X, 'RADIUS  OF  ADDENDUM  CYLINDER  (MM)  RAT2=' ,F14 . 7/ 

&2X, 'RADIUS  OF  DEDENDUM  CYLINDER  (MM)  RDT2=' , F14 . 7/ 

&2X, 'TRANSVERSE  PRESSURE  ANGLE  (RAD.)  PSIT2=' , F14 . 7/ 

&2X, 'HELIX  ANGLE  ON  BASE  CYLINDER(RAD)  LAMDB2=' , F14 . 7/ 

&2X, 'ELASTIC  APPROACH  (MM)  DEL=',F14.7) 

C. 

RETURN 

END 

C . . . 
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c...  ADDITION  OF  TWO  VECTORS 

C. . . 

SUBROUTINE  ADDVEC(VA,VB,VC,DD) 
IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  VA(3) ,VB(3) ,VC(3) 

DO  101  1=1,3 

VA ( I ) =VB (I) +DD*VC (I) 

101  CONTINUE 
RETURN 
END 

C, ,  , 

C. . .  DOT  PRODUCT  OF  TWO  VECTOR 

C. .  . 

SUBROUTINE  DOTVEC(AA, VA,VB) 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  VA(3) ,VB(3) 

AA=O.DO 
DO  102  1=1,3 

AA=AA+VA ( I ) *VB (I) 

102  CONTINUE 
RETURN 
END 

C. 

c'/.  '.  CROSS  PRODUCT  OF  TWO  VECTOR 
C« 

SUBROUTINE  CROVEC(VA,VB,VC) 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  VA(3) ,VB(3) ,VC(3) 

VA(1) =VB { 2 ) *VC (3 ) -VB (3 ) *VC  ( 2 ) 

VA ( 2 ) =VB ( 3 ) * VC ( 1 ) -VB ( 1 ) * VC ( 3 ) 

VA( 3 )  =VB { 1)  *VC (2 ) -VB (2 )  *VC ( 1) 

RETURN 

END 

C. .  . 

c! ! .  PRODUCT  OF  MATRIX  AND  A  VECTOR 

C. .  . 

SUBROUTINE  MAVEC ( VA , MC , VB) 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  MC(3,3) ,VA(3) ,VB(3) 

DO  103  1=1,3 
VA(I)=0.0 
DO  104  J=l,3 

VA(I)=MC(I,J)*VB(J)+VA(I) 

104  CONTINUE 

103  CONTINUE 
RETURN 
END 

C. .  . 

c!  !  PRODUCT  OF  A  VECTOR  AND  A  SCALAR 

C.  .  • 

SUBROUTINE  PDSVEC(VA,VB,T) 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  VA(3) ,VB(3) 

DO  105  1=1,3 
VA(I)=T*VB(I) 

105  CONTINUE 
RETURN 
END 

C. . 

c!  *.  !  STANDARDIZATION  OF  A  VECTOR 
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C.  •  • 

SXJBROUTINE  STDVEC  (VA, VB) 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  VA(3) ,VB(3) 

CC=O.DO 
DO  106  1=1,3 
CC=CC+VB(I)**2 

106  CONTINUE 
CN=DSQRT(CC) 

DO  107  1=1,3 

VA(I)=VB(I)/CN 

107  CONTINUE 
RETURN 
END 

C.  . . 

C.  .  .  INPUT  A  VECTOR  TO  ANOTHER  VECTOR 

C«  •  • 

SUBROUTINE  EQVEC(VA,VB) 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  VA(3),VB{3) 

DO  108  1=1,3 
VA(I)=VB(I) 

108  CONTINUE 
RETURN 
END 

C.  . . 

c! ! !  TRIPLE  PRODUCT  OF  THREE  VECTORS 

C»  •  • 

SUBROUTINE  TRIVEC(AA, VA,VB, VC) 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  VA(3) ,VB(3) ,VC(3) ,V(3) 

CALL  CROVEC(V,VB,VC) 

CALL  DOTVEC(AA,VA,V) 

RETURN 

END 

C  •  •  • 

C...  TRANSFOR  MATRIX 

C. .  • 

SUBROUTINE  TRANSM(AA,BB) 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  AA(3,3) ,BB(3,3) 

AA(1,1)=BB(1,1) 

AA(1,2)=BB(2,1) 

AA(1,3)=BB(3,1) 

AA(2,1)=BB(1,2) 

AA(2,2)=BB(2,2) 

AA{2,3)=BB(3,2) 

AA(3,1)=BB(1,3) 

AA(3,2)=BB(2,3) 

AA{3,3)=BB(3,3) 

RETURN 

END 

C.  . . 

C.  . .  PRODUCT  OF  TWO  MATRICES 

C>  • . 

SUBROUTINE  MAPMA(AA,BB,CC) 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  AA(3,3) ,BB(3,3) ,CC(3,3) 

AA(1,1)=BB(1,1)*CC(1,1)+BB(1,2)*CC(2,1)+BB(1,3)*CC(3,1) 

AA(1,2)=BB(1,1)*CC(1,2)+BB(1,2)*CC(2,2)+BB(1,3)*CC(3,2) 
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AA(1,3)=BB(1,1)*CC(1,3)+BB(1,2)*CC(2,3)+BB(1,3)*CC(3,3) 

AA(2,1)=BB(2,1)*CC(1,1)+BB(2,2)*CC(2,1)+BB(2,3)*CC(3,1) 

AA(2,2)=BB(2,1)*CC(1,2)+BB(2,2)*CC(2,2)+BB(2,3)*CC(3,2) 

AA(2,3)=BB(2,1)*CC(1,3)+BB(2,2)*CC(2,3)+BB(2,3)*CC(3,3) 

AA(3,1)=BB(3,1)*CC(1,1)+BB(3,2)*CC(2,1)+BB(3,3)*CC(3,1) 

AA{3,2)=BB(3,1)*CC(1,2)+BB(3,2)*CC(2,2)+BB(3,3)*CC(3,2) 

AA(3,3)=BB(3,1)*CC(1,3)+BB(3,2)*CC(2,3)+BB(3,3)*CC(3,3) 

RETURN 

END 

C. . 

C. . .  THE  SUBROUTINE  IS  FOR  MISALIGMENT 
SUBROUTINE  HIAL(K,DGAHH,LFH) 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  LFH(3,3) 

S3=^DSIN(DGAMH) 

C3=DCOS(DGAMM) 

C  **  FOR  CROSSING  ANGLE  MISALIGNMENT 
IF(K.EQ.l)  THEM 
LFH{l,l)=--l.DO 
LFH(1,2)*0.D0 
LFH(l,3)=O.DO 
LFH(2,1)=0.D0 
LFH(2,2)=-C3 
LFH{2,3)=-S3 
LFH(3,1)=0.D0 
LFH(3,2)=-S3 
LFH(3,3)*C3 
ENDIF 

C  **  FOR  INTERSECTING  ANGLE  MISALIGNMENT 
IF(K.EQ.2)  THEN 
LFH(1,1)«-C3 
LFH(l,2)=O.DO 
LFH(1,3)=S3 
LFH{2,1)=O.DO 
LFH(2,2)=-1.D0 
LFH(2,3)=0.D0 
LFH(3,1)=S3 
LFH(3,2)=0.D0 
LFH(3,3)=C3 
ENDIF 

C 

RETURN 

END 

C 

C  ....  COMPUTE  THE  PRINCIPAL  CURVATURES  OF  MODIFIED  cutter  SURFACE 

C 

SUBROUTINE  CXJRVT(KK) 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  EFN(3)  ,EHN(3)  ,W1VT2(3)  ,WV12(3)  ,W2VT1(3)  ,KF,KH,KS,KQ 
REAL*8  KM2,KT2,KFF,KHH 

COMMON/A200/  Wl{3) ,W2(3) ,W12(3)  ,VT1(3) ,VT2(3) ,V12(3) 

COMMON  /A210/  EX(3) ,EF(3) ,EH(3) 

COMMON  /A220/  KF,KH 
COMMON  /A230/  ET{3) ,EM(3) 

COMMON  /A300/  ES(3),EQ(3) 

COMMON  /A310/  KFF,KHH 

COMMON  /A380/  KS,KQ 

EFN(1)=  EX(2)*EF(3)-EX(3)*EP(2) 

EFN(2)=>(EX(1)*EF{3)-EX{3)*EF(1)) 

EFN(3)=  EX(1)*EF(2)-EX(2)*EF(1) 
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EHN(1)=  EX(2)*EH{3)-EX(3)*EH(2) 
EHN(2)*-(EX(1)*EH(3)-EX(3)*EH(1)) 

EHN(3)=  EX(1)*EH(2)-EX(2)*EH(1) 

W1VT2(1)*  W1(2)*VT2(3)-W1(3)*VT2(2) 

W1VT2 (2) =- (W1 ( 1) *VT2 ( 3) -W1 (3 ) *VT2 ( 1) ) 

W1VT2(3)=  W1(1)*VT2(2)“W1(2)*VT2(1) 

W2VT1(1)=  W2(2)*VT1(3)-W2(3)*VT1(2) 
W2VT1(2)=-(W2{1)*VT1(3)-W2(3)*VT1(1) ) 

W2VT1(3)=  W2(1)*VT1(2)-W2(2)*VT1(1) 

WV12(1)=  W12(2)*V12(3)-W12(3)*V12(2) 

WV12 ( 2 )  *- (W12 (1) *V12 ( 3 ) -W12 ( 3 )  *V12 ( 1) ) 

WV12 ( 3 ) =  W12 ( 1 ) *V12 ( 2 ) -W12 ( 2 ) *V12 ( 1 ) 

V12F=0 . 0 
V12H=0.0 
WNEF*0.0 
WNEH=0 . 0 
VWN=  0.0 
W1TN=0 . 0 
W2TN=0.0 
DO  1  1=1,3 

V12F=  V12(I)*EF(I)+V12F 
V12H=  V12(I)*EH(I)+V12H 
WNEF=  W12(I)*EFN(I)+WNEF 
WNEH=  W12(I)*EHN(I)+WNEH 
VWN  =  EX(I)*WV12(I)+VWN 
W1TN=  EX(I)*W1VT2(I)+W1TN 

1  W2TN=  EX  (I)  *W2VT1  (I)  +W2TN 

C...  COMPOTE  THE  CURVATURE  OF  THE  GENERATED  SURFACE 
B13=-KF*V12F-WNEF 
B2  3=-KH*V12H-WNEH 

B33=-KF*V12F**2-KH*V12H**2+VWN-W1TN+W2TN 

B11=B13**2/B33 

B12=B13*B23/B33 

B22=B23**2/B33 

Tl=2 . ODO  0  *B13  *B2  3 

T2=B23**2-B13**2-(KF-KH) *B33 

SIG1F=0 . 5DOO*DATAN (T1/T2 ) 

C. ..  .PRINCIPAL  CURVATURES  OF  THE  GENERATED  SURFACE 
IF(DABS(SIG1F) .LE.O.lD-5)  THEN 

KQ=0.5D0*((KF+KH)+(B13**2+B23**2)/B33 
&  +(B23**2-B13**2-(KF-KH) *B33)  / (B33*DCOS (SIGIF) ) ) 

KS=KQ-(B23**2-B13**2-(KF-KH)*B33)/(B33*DC0S(SIG1F)  ) 
ELSE 

KQ=0 . 50D00* (KF+KH) +0 . 5D00* (B13**2+B23**2) /B33 
&  +B13*B23/ (B33*DSIN(2 . ODOO* (+SIG1F) ) ) 

KS=  KQ-2.0D00*B13*B23/ (B33*DSIN(2. ODOO* (+SIG1F) ) ) 
ENDIF 

SIGSF=-SIG1F 

C.  .  .  .PRINCIPW:.  DIRECTIONS  OF  THE  GENERATED  SURFACE 
DO  2  1=1,3 

EQ ( I ) =  DCOS ( SIGSF) *EH ( I ) -DSIN { SIGSF) *EF ( I ) 

2  ES(I)=  DSIN(SIGSF)*EH(I)+DCOS(SIGSF)*EF(I) 

IF  (KK  .LT.  2)  GO  TO  100 

Q3=ES(1)*ET(1)+ES(2)*ET(2)+ES(3)*ET(3) 

Q4=EQ(1)*ET(1)+EQ(2)*ET(2)+EQ(3)*ET(3) 

IF(Q3.EQ.Q4)  THEN 
Q34=DATAN(Q4) 

ELSE 

Q34=DATAN2(Q4,Q3) 

ENDIF 
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KT2=KS*DCOS(Q34) **2+KQ*DSIN(Q34) **2 
KM2=KS*DSIN(Q34) **2+KQ*DC0S (Q34)**2 
GO  TO  120 
100  KFF=KS 
KHH=KQ 
120  RETURN 
END 
C 

C  ....  COMPUTE  THE  DIMENSIONS  AND  DIRECTIONS  OF  THE  CONTACT  ELLIPSE 
C 

SUBROUTINE  ELLIP 
IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  KSIG1,KSIG2,KFF,KHH,KS,KQ 
COMMON  /A300/  ES(3),EQ(3) 

COMMON  /A380/  KS,KQ 
COMMON  /A310/  KFF,KHH 
COMMON  /A340/  EFF(3) ,EHH(3) 

COMMON  /A360/  A,B,SI(3) ,FI(3) 

DEL=0.007D0 

PI®DACOS(-1.DO) 

S1=ES ( 1) *EFF (1) +ES ( 2 ) *EFF ( 2 ) +ES (3 ) *EFF (3 ) 
S2=EQ{1)*EFF(1)+EQ(2)*EFF(2)+EQ(3)*EFF(3) 

SIGSF=DATAN2 (S2 , SIJ 

C.... COMPUTE  THE  DIMENSIONS  OF  THE  CONTACT  ELLIPSE  (A  &  B) 
KSIG1=KFF+KHH 
KSIG2=KS+KQ 
G1=KFF-KHH 
G2=KS-KQ 

A=  (KSIG1-KSIG2-  (Gl**2-2 .  DO*G1*G2*DCOS  (2 .  D0*SIGSF)  +G2**2 )  **0 . 5) 
&  /4.D0 

B=(KSIG1-KSIG2+(G1**2-2.DO*G1*G2*DCOS(2.DO*SIGSF)+G2**2)  **0.5) 
&  /4.D0 

A= (DEL/ABS (A) ) **0 . 5 
B=(DEL/ABS(B) )**0.5 

C.... COMPUTE  THE  ANGLE  (ALFl)  BETWEEN  AXES  OF  ELLIPSE  &  PRINCIPLE 
S1=G2  *DSIN ( 2 . D0*SIGSF) 

S2=G1-G2*DC0S (2 . D0*SIGSF) 

ALF1=0.5D0*DATAN2 (S1,S2) 

S3=DSQRT(G1**2-2.DO*G1*G2*DCOS(2.DO*SIGSF)+G2**2) 

SS2=S1/S3 

SC2=S2/S3 

ALF1=DATAN(SS2/ (1.D0+SC2) ) 

ALF1=DABS (ALFl) 

C..  AXES  OF  THE  CONTACT  ELLIPSE 
DO  100  1=1,3 

SI ( I ) =DSIN (ALFl ) *  EFF ( I ) +DCOS ( ALFl ) *EHH ( I ) 

FI (I) =DCOS (ALFl) *EFF(I) “DSIN(ALFl) *EHH(I) 

100  CONTINUE 
RETURN 
END 
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